Fault location using measurements of current and voltage from one end of a line

ABSTRACT

A method to locate a fault from one end of a section of a power line utilizing measurements of current, voltage and angles between the phases at a first end of said section. Symmetrical components of currents are calculated for the current and voltage measurement at the first end. A value of impedance is calculated for an extra link between the terminals with the impedance for the positive sequence. A compensation is determined for the shunt capacitance. The zero-sequence current is determined from the healthy line of a section of parallel power lines. A distance to a fault is calculated for the parallel line section. The distance to the fault from the first end is calculated. The fault is located utilizing the calculate distances.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority to Swedish patent application 0201911-5 filed 20 Jun. 2002 and is the national phase under 35 U.S.C. § 371 filed 18 Jun. 2003.

TECHNICAL FIELD

The present invention is concerned with a fault location technique for a section of a power transmission line utilizing measurements of current and voltage made at terminals located at one end of the section of the power line.

BACKGROUND ART

Several methods and approaches for fault location in high voltage power transmission systems, and power distribution systems, have been developed and employed. One approach has been to use voltage/current transducers located at terminals located at each of two ends of a section of the power line to be monitored. Inductive current transformers are used to provide a measurement of instantaneous current in a transmission line.

U.S. Pat. No. 4,559,491 entitled Method and device for locating a fault point on a three-phase transmission line, describes a one-end fault location (FL) algorithm. High accuracy of fault location using a fault locator device at one end of a line is achieved by taking into account the actual distribution of a fault current in the transmission network. This algorithm has been successfully implemented into a product in 1982 and is in operation with single and parallel transmission lines in many countries around the world. However, for certain fault conditions it is difficult to obtain accurate pre-fault quantities, such as pre-fault currents, in order to calculate an estimate for voltage drop across the fault path. Also, a disadvantage of using phase voltages and currents and zero sequence components of currents is that it is relatively difficult using these values to compensate for shunt disadvantage of using phase voltages and currents and zero sequence components of currents is that it is relatively difficult using these values to compensate for shunt capacitance effects. In addition, the fault locator method described is not suitable for single and parallel line sections which have an extra link across the ends of the sections.

SUMMARY OF THE INVENTION

The aim of the present invention is to remedy one or more of the above mentioned problems.

This is obtained by a method characterised by claim 1. Specific features of the present invention are characterised by the appending claims.

In one aspect of the invention, a method comprising a new formulation of a one-end fault locator algorithm has been proposed. The uniform description of the transmission network in terms of symmetrical components as well as the generalized models of fault loops and faults have been applied. The resulting advantages include the algorithm can be used for locating faults in typical single and parallel transmission lines, and, in addition, fault location may also be carried out for both single and parallel lines with an extra link between the line ends. Another advantage is that a procedure for calculation of a distance to fault is in the form of a compact quadratic equation with the coefficients dependent on a fault type, acquired measurements and impedance data for the transmission network. Another advantage of the invention is that optimal estimation of the voltage drop across a fault path is applied, which has the result that the pre-fault currents in case of single phase-to-ground faults and phase-to-phase faults are no longer required.

In an embodiment, compensation for shunt capacitances is facilitated by means of the use of the notation of symmetrical components. The distributed long line model of the line has been applied for that. The compensation is performed individually for all the sequences. The currents for particular sequences are compensated for the shunt currents and then the fault loop compensated current is composed. In another embodiment improved accuracy has been obtained by means of an option to measure the source impedance at the remote end instead of using a representative value. The source impedance measured at the remote end may be considered as sent to the fault locator by using a simple communication means.

In another embodiment, a method for one end fault location for parallel lines to locate single phase-to-ground faults is described under a plurality of conditions. In another further embodiment a method is described for one end fault location with standard availability of the measured signals for ground faults including both single phase-to-ground faults and phase-to-phase-to-ground faults.

In another aspect of the invention, a fault locator device for carrying out the method of the invention is characterised by claim 18. Specific features of the inventive fault locator device are characterised by the appending claims.

In another aspect of the invention a computer program is described for carrying out the method according to the invention. In another aspect of the invention a computer program product comprising a computer program for carrying out the method of the invention is described. In another, further aspect of the invention a graphical user interface is described for displaying a distance to a fault from one end of a section of a power line.

BRIEF DESCRIPTION OF THE DRAWINGS

A more complete understanding of the method and system of the present invention may be had by reference to the following detailed description when taken in conjunction with the accompanying drawings wherein:

FIG. 1 shows in a single schematic diagram a method of fault location in power transmission and/or distribution systems for parallel lines and single lines according to an embodiment of the invention;

FIG. 2 a shows a schematic circuit diagram for a parallel transmission network for the positive sequence component in which the fault loop is marked for the case of a fault locator installed at the terminal AA. FIGS. 2 b, 2 c show corresponding diagrams for the negative sequence and zero sequence components, respectively;

FIG. 3 a is a schematic block diagram for obtaining and calculating the phasors of the symmetrical components of voltages and currents used for composing the fault loop voltage. FIG. 3 b shows a corresponding diagram for composing the fault loop current;

FIG. 4 shows a circuit diagram for determining the fault current distribution factor for the positive sequence of a single line, in which diagram quantities for the negative sequence are shown indicated in brackets;

FIG. 5 shows a circuit diagram corresponding to FIG. 4 for single lines for determining the fault current distribution factor for the positive sequence of parallel lines, in which quantities for the negative sequence are also shown indicated in brackets;

FIG. 6 shows a schematic diagram for an embodiment of the invention in which source impedance measured at a remote end B may be communicated to a fault locator at the first end A;

FIG. 7 is a circuit diagram of an embodiment in which the shunt capacitances are taken into account, and shows a positive sequence circuit diagram during a first iteration;

FIG. 8 shows a negative sequence circuit diagram for taking the shunt capacitances effect into account during a first iteration;

FIG. 9 shows a zero sequence circuit diagram for taking the shunt capacitances effect into account during a first iteration;

FIG. 10 shows a flowchart for a method for locating a fault in a single line according to an embodiment of the invention;

FIG. 11 shows a flowchart for a method for locating a fault in parallel lines according to an embodiment of the invention;

FIG. 12 and FIGS. 13 a, 13 b, 14, 15 a and 15 b show schematic diagrams of possible fault-types (phase-to-phase, phase-to-ground and so on) with respect to derivation of coefficients for Table 2 in Appendix A2. FIG. 12 shows fault types from a-g, and FIGS. 13 a, 13 b faults between phases a-b. FIG. 14 shows an a-b-g fault. FIGS. 15 a and 15 b show symmetrical faults a-b-c and a-b-c-g respectively;

FIGS. 16 and 17 show schematic diagrams for the derivation of the complex coefficients in the fault current distribution factors for the positive (negative) sequence included in Table 3. FIG. 16 shows the case of a single line with an extra link between the substations. FIG. 17 shows the case of parallel lines with an extra link between the substations;

FIG. 18 shows a fault locator device and system according to an embodiment of the invention;

FIG. 19 shows a flowchart for a method for locating a single phase-to-ground fault in parallel lines in the case of measurements from the healthy line being unavailable, according to an embodiment of the invention;

FIG. 20 shows a schematic diagram for a method of fault location for parallel lines with different modes of the healthy parallel operation;

FIG. 21 a, shows a schematic equivalent circuit diagram for a parallel network for the incremental positive or the negative sequence.

FIG. 21 b shows an equivalent circuit diagram for the zero sequence while both parallel lines are in operation.

FIG. 21 c shows the equivalent circuit diagram for the zero sequence with the healthy parallel line switched off and grounded;

FIG. 22 shows a flowchart for a method for locating phase-to-phase and phase-to-ground faults in parallel lines in the case of providing the zero sequence currents from the healthy parallel line according to another embodiment of the invention;

FIG. 23 shows a schematic diagram fault location for parallel lines with standard availability of measurements according to another embodiment of the invention;

FIGS. 24 a,b,c show the equivalent circuit diagrams of parallel lines for positive, negative and zero sequence currents respectively.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

FIG. 1 presents a schematic diagram for one-end fault location applied for parallel lines and for a power transmission or distribution system with a single line. A fault locator 1 is positioned at one end 2 of a single line AA-BA 3 or parallel lines AA-BA, AB-BB, 4. A fault F is shown at FA with a corresponding fault resistance, 5, denoted as R_(F). A value for distance to the fault d from one end 2 determined and provided by the fault locator 1 is indicated with the reference number 7. Components such as parallel line AB-BB and quantities such as a parallel line value zero sequence current I_(AB0) shown with dotted lines are excluded when considering a single line case.

The fault locator 1 positioned at the first end 2, or ‘A’ end, is supplied with the following input signals:

-   -   three-phase voltages (V_(AA)) of the faulted line,     -   three-phase voltages (I_(AA)) of the faulted line,     -   zero sequence current (I_(AB0)) from the healthy parallel line         (zero sequence is not present when the single line only is         considered).

FIGS. 2 a,b,c show circuit diagrams of a parallel transmission network for the positive 2 a, negative 2 b, and zero sequence 2 c components. Fault loops for the sequence components 21 a, 21 b, 21 c are shown for the case of a fault locator installed at the terminal AA. An extra link 25 between the terminals A, B is shown. A generalized model of the fault loop considered for different fault types is stated as: V _(AA) _(—) _(p) −dZ _(1LA) I _(AA) _(—) _(p) −R _(F)(a _(F1) I _(F1) +a _(F2) I _(F2) +a _(F0) I _(F0))=0   (1) where: d—unknown and sought distance to fault, Z_(1LA)—positive sequence impedance of the faulted line, V_(AA) _(—) _(p), I_(AA) _(—) _(p)—fault loop voltage and current composed according to the fault type, R_(F)—fault resistance, I_(Fi)—sequence components of the total fault current (i=0, i=1, i=2), a_(Fi)—weighting coefficients (TABLE 2).

Fault loop voltage and current can be expressed in terms of the symmetrical components of measured voltages/currents: V _(AA) _(—) _(p) =a ₁ V _(AA1) +a ₂ V _(AA2) +a ₀ V _(AA0)   (2)

$\begin{matrix} {{\underset{\_}{I}}_{{AA}{\_ p}} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{AA2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{AA0}} + {{\underset{\_}{a}}_{0m}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{La}}}{\underset{\_}{I}}_{AB0}}}} & (3) \end{matrix}$ where: AA, AB—subscripts used for indicating measurements acquired from the faulted line (AA) and from the healthy line (AB), respectively, a₀, a₁, a₂—coefficients which are gathered in TABLE 1 (the Tables are arranged below at the end of the description of embodiments and derivation of these coefficients is shown in APPENDIX A1, also attached). Z_(0LA), Z_(0m)—impedance of the faulted line and mutual coupling between the lines for the zero sequence, respectively, a_(0m)=a₀—for parallel lines, a_(0m)=0—for single lines.

The phasors of symmetrical components of voltages, positive: V_(AA1), negative: V_(AA2) and zero sequence: V_(AA0) as well as the phasors of symmetrical components of currents, positive: I_(AA1), negative: I_(AA2), zero sequence from the faulted line: I_(AA0) and zero sequence from the healthy line: I_(AB0) are calculated from the acquired measurements as shown schematically in schematic block diagrams FIGS. 3 a and 3 b.

FIG. 3 a shows an input of instantaneous phase voltages 30 a, filtering stage 33 a, phasors of phase voltages 31 a, calculation of phasors of symmetrical components 33 b and phasors of symmetrical components of voltages output at 32 a. It may be seen from FIG. 3 a that acquired phase voltage measurements are subjected to a filter, then calculations are made to find the symmetrical components of the fault loop voltage. FIG. 3 b shows correspondingly stages used to find the symmetrical components of the fault loop current. FIG. 3 b shows instantaneous phase currents and instantaneous zero sequence current from the healthy line 30 b, filtering 33 b, phasors of phase currents and phasor of zero sequence current from the healthy line 31 b, calculation 34 b and phasors of symmetrical components of currents output at 32 b.

Fault loop signals may be composed according to formulae (2)-(3) and TABLE 1, which is the alternative to the classic approach (TABLE 1A, fault loop voltage (V_(AA) _(—) _(p)) and current (I_(AA) _(—) _(p)), which was used in the fault locator from [1-2].

Voltage drop across a fault path resistance, the third term in (1), can be expressed in terms of the current distribution factors and local measurements of currents which results in:

$\begin{matrix} {{{\underset{\_}{V}}_{AA\_ p} - {d{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ p}} - {R_{F}\left( {{{\underset{\_}{a}}_{F1}\frac{\Delta\;{\underset{\_}{I}}_{AA1}}{{\underset{\_}{k}}_{F1}}} + {{\underset{\_}{a}}_{F2}\frac{{\underset{\_}{I}}_{AA2}}{{\underset{\_}{k}}_{F2}}} + {{\underset{\_}{a}}_{F0}\frac{{\underset{\_}{I}}_{AA0}}{{\underset{\_}{k}}_{F0}}}} \right)}} = 0} & (4) \end{matrix}$

Formula (4) has been obtained from the following relations between the symmetrical components of a total fault current and measured currents:

$\begin{matrix} {{{\underset{\_}{I}}_{F1} = \frac{\Delta\;{\underset{\_}{I}}_{AA1}}{{\underset{\_}{k}}_{F1}}},\mspace{14mu}{{\underset{\_}{I}}_{F2} = \frac{{\underset{\_}{I}}_{AA2}}{{\underset{\_}{k}}_{F2}}},\mspace{14mu}{{\underset{\_}{I}}_{F0} = \frac{{\underset{\_}{I}}_{AA0}}{{\underset{\_}{k}}_{F0}}}} & (5) \end{matrix}$ where: I_(F1); I_(F2); I_(F0)—symmetrical components of a total fault current, k_(F1), k_(F2), k_(F0)—fault current distribution factors for particular sequence quantities, ΔI_(AA1)=I_(AA1)−I_(AA1pre); I_(AA2); I_(AA0)—symmetrical components of currents measured in the line A at the station A (subscript AA); note that in case of the positive sequence the incremental quantity (post-fault current minus pre-fault current) is used.

Voltage drop across the fault path, as shown in the third term in equation (1), is expressed using sequence components of the total fault current. The weighting coefficients a_(F0), a_(F1), a_(F2), can accordingly be determined by taking the boundary conditions for particular fault type. See TABLE 2, Alternative sets of the weighting coefficients for determining a voltage drop across the fault path resistance. Examples of derivation of these coefficients are contained in APPENDIX A2.

There is some freedom for setting the weighting coefficients. It is proposed to utilize this freedom firstly for avoiding zero sequence quantities, since the zero sequence impedance of a line may be considered as an unreliable parameter. This can be accomplished by setting a_(F0)=0 as shown in TABLE 2.

Secondly, the freedom in establishing the weighting coefficients can be utilized for determining the preference for using particular quantities. The negative sequence (TABLE 2, set I) or the positive sequence (TABLE 2, set II) can be preferred as well as possibly both types of the quantities (TABLE 2, set III) can be used for determining the voltage drop across the fault path.

The set I is recommended for further use, thus avoiding the positive sequence, and thus avoiding the pre-fault positive sequence current, for the largest number of faults. Avoiding the pre-fault positive sequence current is highly desirable since sometimes the pre-fault currents—due to certain reasons—can not be recorded or registered, but may be contaminated by one or more the symptoms of the occurring fault. Moreover, the accuracy of recording the pre-fault currents, which are basically lower than the post-fault currents, is not very great. This is so since the A/D converters operate with less accuracy in the low range.

Fault current distribution factors depend on the configuration of the transmission network, FIGS. 4, 5, and impedance parameters. Basically, all impedances for the positive and for the negative sequence are equal to each other and thus one obtains:

$\begin{matrix} {{\underset{\_}{k}}_{F1} = {{\underset{\_}{k}}_{F2} = \frac{{{\underset{\_}{K}}_{1}d} + {\underset{\_}{L}}_{1}}{{\underset{\_}{M}}_{1}}}} & (6) \end{matrix}$ Coefficients in a fault current distribution factor (6) for a single (FIG. 4) and for parallel lines (FIG. 5) are gathered in TABLE 3, Coefficients for determining a fault current distribution factor, (note that derivation of the coefficients is shown in APPENDIX A3).

FIG. 4 shows a circuit diagram of a single line for determining the fault current distribution factor for the positive sequence currents and with the negative sequence currents such as shown in brackets. Similarly FIG. 5 shows a circuit diagram of parallel lines for determining the fault current distribution factor with positive sequence currents wherein the negative sequence currents are shown in brackets.

In FIG. 4 the extra link 45 between the terminals A, B having impedance for the positive sequence equal to: Z_(1AB) can be considered as existing (Z_(1AB)≠∝) or as not present (Z_(1AB)→∝). In FIG. 5 the extra link 55 between the terminals A, B having impedance for the positive sequence equal to: Z_(1AB) can be considered as existing

$\left( {{\underset{\_}{Z}}_{{{1{LB}}\&}{AB}} = \frac{{\underset{\_}{Z}}_{1{LB}}{\underset{\_}{Z}}_{1{AB}}}{{\underset{\_}{Z}}_{1{LB}} + {\underset{\_}{Z}}_{1{AB}}}} \right)$ or as not present (Z_(1LB&AB)=Z_(1LB)).

Substituting (6) into (4) and adjusting a_(F0)=0 (as in TABLE 2) results in:

$\begin{matrix} {{{\underset{\_}{V}}_{AA\_ p} - {d{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ p}} - {\frac{R_{F}{\underset{\_}{M}}_{1}}{{{\underset{\_}{K}}_{1}d} + {\underset{\_}{L}}_{1}}\left( {{{\underset{\_}{a}}_{F1}\Delta\;{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{AA2}}} \right)}} = 0} & (7) \end{matrix}$

After multiplying both sides of (7) by:

$\frac{{{\underset{\_}{K}}_{1}d} + {\underset{\_}{L}}_{1}}{{\underset{\_}{I}}_{AA\_ p}}$ and some rearrangements, the quadratic formula with two unknowns, d-[p.u.] sought fault distance from A, R_(F)—fault resistance, is obtained:

$\begin{matrix} {{{{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{1L}d^{2}} + {\left( {{{\underset{\_}{L}}_{1}{\underset{\_}{Z}}_{1L}} - {{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{AA\_ p}}} \right)d} - {{\underset{\_}{L}}_{1}{\underset{\_}{Z}}_{AA\_ p}} + {R_{F}{\underset{\_}{M}}_{1}\frac{\left( {{{\underset{\_}{a}}_{F1}\Delta\;{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{AA2}}} \right)}{{\underset{\_}{I}}_{AA\_ p}}}} = 0.} & (8) \end{matrix}$ where:

${\underset{\_}{Z}}_{AA\_ p} = {\frac{{\underset{\_}{V}}_{AA\_ p}}{{\underset{\_}{I}}_{AA\_ p}} - {{calculated}\mspace{14mu}{fault}\mspace{14mu}{loop}\mspace{14mu}{{impedance}.}}}$

Writing formula (8) in more compact form results in: A ₂ d ² +A ₁ d+A ₀ +A ₀₀ R _(F)=0  (8a) where: A ₂ =A ₂ _(—) _(Re) +jA ₂ _(—) _(Im) =K ₁ Z _(1LA) A ₁ =A ₁ _(—) _(Re) +jA ₁ _(—) _(Im) =L ₁ Z _(1LA) −K ₁ Z _(AA) _(—) _(p) A ₀ =A ₀ _(—) _(Re) +jA ₀ _(—) _(Im) =−L ₁ Z _(AA) _(—) _(p)

$A_{00{\_ Re}} = {{j\; A_{oo\_ Im}} = \frac{{\underset{\_}{M}}_{1}\left( {{{\underset{\_}{a}}_{F1}\Delta\;{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{AA2}}} \right)}{{\underset{\_}{I}}_{AA\_ p}}}$

${\underset{\_}{Z}}_{AA\_ p} = \frac{{\underset{\_}{V}}_{AA\_ p}}{{\underset{\_}{I}}_{AA\_ p}}$

-   -   calculated fault loop impedance K₁, L₁, M₁—coefficients gathered         in TABLE 3.

Formula (8a) can be written separately for real and imaginary parts: A ₂ _(—) _(Re) d ² +A ₁ _(—) _(Re) d+A ₀ _(—) _(Re) +A ₀₀ _(—) _(Re) R _(F)=0  (8b) A ₂ _(—) _(Im) d ² +A ₁ _(—) _(Im) d+A ₀ _(—) _(Im) +A ₀₀ _(—) _(Im) R _(F)=0  (8c)

Combining (8b) and (8c) in such the way that fault resistance is eliminated [that is, equation (8b) is multiplied by A₀₀ _(—) _(Im) and equation (8c) by A₀₀ _(—) _(Re) and then subtracting them] yields the quadratic formula for a sought fault distance: B ₂ d ² +B ₁ d+B ₀=0  (9) where: B ₂ =A ₂ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₂ _(—) _(Im) A ₀₀ _(—) _(Re) B ₁ =A ₁ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₁ _(—) _(Im) A ₀₀ _(—) _(Re) B ₀ =A ₀ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₀ _(—) _(Im) A ₀₀ _(—) _(Re) Equation (9) has two roots (d₁, d₂) for the distance to fault:

$\begin{matrix} {{d_{1} = \frac{{- B_{1}} - \sqrt{B_{1}^{2} - {4B_{2}B_{0}}}}{2B_{2}}}{d_{2} = \frac{{- B_{1}} + \sqrt{B_{1}^{2} - {4B_{2}B_{0}}}}{2B_{2}}}} & (10) \end{matrix}$ The root which fulfils the condition (0≦d≦1) is selected as the solution for the distance to fault.

In another embodiment of the invention, the method of fault location is carried out by using a measurement of source impedance at the second end remote from the fault locator 1, instead of a representative value for the source impedance at the remote end, and by communicating that measured value to the local end using a communication means. Coefficients from (9) are determined with the local measurements and the impedance data for the transmission line, the extra link between the line terminals and the equivalent systems at the line terminals. Impedance of the equivalent system at the local substation (Z_(1sA)) can be traced on-line with the local measurements. In contrast, the remote system impedance (Z_(1sB)) is not measurable locally from A. Thus, the “representative” value of this impedance may be provided for the algorithm [1-2].

The alternative solution for the single line case is shown in FIG. 6, which shows a fault locator 1 at the first end 2 near to system A, and another device 10 located at the remote end 3 close to system B, indicated as RD. A communication signal 9 is shown being sent from the 10 at the remote end to the fault locator 1 at the local end.

The remote source impedance (Z_(1sB)) is measured by the remote device RD, 10, which may be another fault locator or any suitable device such as a digital relay or digital fault recorder, in the remote substation and the measurement 9 is sent via a communication channel 60. Synchronization of measurements at the line terminals is not required. The source impedance is calculated from the known relation between the incremental positive voltage (ΔV_(B1)) and the incremental positive sequence current (ΔI_(B1)) [3-4]:

$\begin{matrix} {{\underset{\_}{Z}}_{1{sB}} = {- \frac{\Delta\;{\underset{\_}{V}}_{B1}}{\Delta\;{\underset{\_}{I}}_{B1}}}} & (11) \end{matrix}$

Similarly, the fault locator 1 calculates the local source impedance

$\begin{matrix} {{\underset{\_}{Z}}_{1{sA}} = {- \frac{\Delta\;{\underset{\_}{V}}_{A1}}{\Delta\;{\underset{\_}{I}}_{A1}}}} & \left( {11a} \right) \end{matrix}$

In another and preferred embodiment of the invention, compensation is carried out for the shunt capacitance of the line. Compensation for shunt capacitances effect can be accomplished by taking into account the lumped model (only the longitudinal R-X_(L) parameters are taken into account) or the distributed long transmission line model. The distributed long line model [5] as providing higher accuracy of fault location, has been considered here.

The compensation for the single line is presented further. This means that when composing fault loop current (3) the term reflecting the mutual coupling effect disappears (a_(0m)=0). Moreover, the single subscript (A instead of AA) is used.

Fault location procedure with compensating for shunt capacitances of a transmission line requires the following additional input data:

C_(1L)—shunt capacitance of a whole line for the positive and the negative sequences (parameters of a line for the positive and the negative sequences are identical and thus: C_(2L)=C_(1L))

C_(0L)—shunt capacitance of a whole line for the zero sequence,

l—total line length (km), used for expressing impedances/capacitances of the line per km length.

The compensation of shunt capacitances may be introduced while determining the voltage drop across the faulted line segment—the second term in the generalized fault loop model (1). This requires compensating the components of the calculated currents for particular sequences. Thus, the original measured currents: I_(A1), I_(A2), I_(A0) have to be replaced by the currents after the introduced compensation: I_(A1) _(—) _(comp), I_(A2) _(—) _(comp), I_(A0) _(—) _(comp). At the same time the original fault loop voltage, the first term in the model (1), is taken for a distance to fault calculation. As concerns determining the voltage drop across the fault resistance, the third term in (1), it is assumed here, which is a standard practice, that the effect of line capacitances at the location of the fault (point F), may be neglected. This is justified as the impedance of the capacitive branch at that location is much greater than the fault resistance. This means that the voltage drop across the fault resistance is determined without taking into account the shunt capacitances.

Calculating a distance to fault the following impedances (defined below) are taken as:

Z_(1L) ^(long)—positive sequence impedance of a line with taking into account the distributed long line model,

Z_(0L) ^(long)—as above, but for the zero sequence.

The compensation procedure requires iterative calculations, performed until the convergence is achieved (i.e. iterated until the location estimate ceases to differ from the previous estimate). However, the studies conducted revealed that results of acceptable accuracy may be obtained using a fixed number of iterations, for example, 2-3 iterations. The calculated distance to a fault from a particular (say, present iteration) is utilized for determining the shunt current in the next iteration. The determined shunt current is then deduced from the measured current. A distance to fault, calculated without considering the shunt effect (10), is taken as the starting value for the first iteration.

A way of conducting the first iteration of the compensation is shown in FIGS. 7, 8, 9 for the positive sequence, negative sequence and zero sequence respectively with taking into account the shunt capacitances effect.

As a result of performing the first iteration for the positive sequence (FIG. 7) the compensated current I_(A1) _(—) _(comp) _(—) ₁ is calculated and the last index in the subscript denotes the first iteration. The calculation is based on deducing the shunt current from the calculated positive sequence current I_(A1) calculated from the measured phase currents—FIG. 2: I _(A1) _(—) _(comp) _(—) ₁ =I _(A1)−0.5d _(v) lB′ _(1L) A _(tanh 1) V _(A1)   (12) where:

d_(v)—distance to fault calculated without taking into account the shunt capacitance effect (10),

l—total line length (km)

${\underset{\_}{A}}_{tanh1} = \frac{\tanh\left( {\sqrt{0.5{\underset{\_}{Z}}_{1L}^{\prime}{\underset{\_}{B}}_{1L}^{\prime}}d_{v}l} \right)}{\sqrt{0.5{\underset{\_}{Z}}_{1L}^{\prime}{\underset{\_}{B}}_{1L}^{\prime}}d_{v}l}$

${\underset{\_}{B}}_{1L}^{\prime} = \frac{j\;\omega\; C_{1L}}{l}$

-   -   positive sequence admittance (capacitive) of a line per km         length (S/km)

${\underset{\_}{Z}}_{1L}^{\prime} = \frac{{\underset{\_}{Z}}_{1L}}{l}$

-   -   positive sequence impedance of a line per km length (Ω/km)

Positive sequence impedance of a faulted line segment, between points A and F, without taking into account the shunt capacitances effect and considering the lumped line model equals: d_(v)lZ′_(1L)  (13) while for the considered here distributed long line model: d_(v)lZ′_(1L)A_(sinh 1)  (14) where:

${\underset{\_}{A}}_{sinh1} = \frac{\sinh\left( {\sqrt{{\underset{\_}{Z}}_{1L}^{\prime}{\underset{\_}{B}}_{1L}^{\prime}}d_{v}l} \right)}{\sqrt{{\underset{\_}{Z}}_{1L}^{\prime}{\underset{\_}{B}}_{1L}}d_{v}l}$

Thus, the positive sequence impedance of a line with taking into account the distributed long line model (Z_(1L) ^(long)) equals: Z_(1L) ^(long)=A_(sinh 1)Z_(1L)  (15)

As a result of performing the first iteration for the negative sequence (FIG. 8) the compensated current I_(A2) _(—) _(comp) _(—) ₁ is calculated and the last index in the subscript denotes the first iteration. This is based on deducing the shunt current from the calculated negative sequence current I_(A2), calculated from the measured phase currents—FIG. 2): I _(A2) _(—) _(comp) _(—) ₁ =I _(A2)−0.5d _(v) lB′ _(2L) A _(tanh 2) V _(A2)  (16) where, taking into account that the line parameters for the positive and for the negative sequences are identical (C_(2L)=C_(1L), Z_(2L)=Z_(1L)): A_(tanh 2)=A_(tanh 1) B′_(2L)=B′_(1L)

As a result of performing the first iteration for the zero sequence (FIG. 9) the compensated current I_(A0) _(—) _(comp) _(—) ₁ is calculated, last index in the subscript denotes the first iteration. This calculation is based on deducing the shunt current from the calculated zero sequence current I_(A0), calculated from the measured phase currents—FIG. 2: I _(A0) _(—) _(comp) _(—) ₁ =I _(A0)−0.5d _(v) lB′ _(0L) A _(tanh 0) V _(A0)  (17) where:

${\underset{\_}{A}}_{tanh0} = \frac{\tanh\left( {\sqrt{0.5{\underset{\_}{Z}}_{0L}^{\prime}{\underset{\_}{B}}_{0L}^{\prime}}d_{v}l} \right)}{\sqrt{0.5{\underset{\_}{Z}}_{0L}^{\prime}{\underset{\_}{B}}_{0L}^{\prime}}d_{v}l}$

${\underset{\_}{B}}_{0L}^{\prime} = \frac{j\;\omega\; C_{0L}}{l}$

-   -   zero sequence admittance (capacitive) of a line per km length         (S/km)

${\underset{\_}{Z}}_{0L}^{\prime} = \frac{{\underset{\_}{Z}}_{0L}}{l}$

-   -   zero sequence impedance of a line per km length (Ω/km)

Zero sequence impedance of a faulted line segment, that is between points A and F, without taking into account the shunt capacitances effect and considering the lumped line model equals: d_(v)lZ′_(0L)  (18) while for the considered here distributed long line model: d_(v)lZ′_(0L)A_(sinh 0)  (19) where:

${\underset{\_}{A}}_{sinh0} = \frac{\sinh\left( {\sqrt{{\underset{\_}{Z}}_{0L}^{\prime}{\underset{\_}{B}}_{0L}^{\prime}}d_{v}l} \right)}{\sqrt{{\underset{\_}{Z}}_{0L}^{\prime}{\underset{\_}{B}}_{0L}^{\prime}}d_{v}l}$

Thus, the zero sequence impedance of a line with taking into account the distributed long line model (Z_(0L) ^(long)) equals: Z_(0L) ^(long)=A_(sinh 0)Z_(0L)  (20)

The quadratic complex formula (8) with two unknowns (d_(comp) _(—) ₁[p.u.]—sought fault distance, R_(F)—fault resistance) after introducing the compensation (first iteration) takes the following form:

$\begin{matrix} {{{{{\underset{\_}{K}}_{1}{{\underset{\_}{Z}}_{1L}^{long}\left( d_{{comp\_}1} \right)}^{2}} + {\left( {{{\underset{\_}{L}}_{1}{\underset{\_}{Z}}_{1L}^{long}} - {{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{A\_ p}^{{comp\_}1}}} \right)d_{{comp\_}1}} - {{\underset{\_}{L}}_{1}{\underset{\_}{Z}}_{A\_ p}^{{comp\_}1}} + {R_{F}{\underset{\_}{M}}_{1}\frac{\left( {{{\underset{\_}{a}}_{F1}\Delta\;{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{AA2}}} \right)}{{\underset{\_}{I}}_{A\_ p}^{{comp\_}1}}}} = 0}{{where}\text{:}}{{{\underset{\_}{Z}}_{1L}^{long} = {{\underset{\_}{A}}_{sinh1}{\underset{\_}{Z}}_{1L}}},{{\underset{\_}{A}}_{sinh1} = \frac{\sinh\left( {\sqrt{{\underset{\_}{Z}}_{1L}^{\prime}{\underset{\_}{B}}_{1L}^{\prime}}d_{v}l} \right)}{\sqrt{{\underset{\_}{Z}}_{1L}^{\prime}{\underset{\_}{B}}_{1L}^{\prime}}d_{v}l}},}} & (21) \end{matrix}$

${\underset{\_}{Z}}_{A\_ p}^{{comp\_}1} = \frac{{\underset{\_}{V}}_{A\_ p}}{{\underset{\_}{I}}_{A\_ p}^{{comp\_}1}}$

-   -   fault loop impedance calculated with:         V_(A) _(—) _(p)—original fault loop voltage (2),         I_(A) _(—) _(p) ^(comp) ^(—) ¹=a₁I_(A1) _(—) _(comp) _(—)         ₁+a₂I_(A2) _(—) _(comp) _(—) ₁+a₀I_(A0) _(—) _(comp) _(—)         ₁—fault loop current composed of the positive (12),         negative (16) and zero (17) sequence currents obtained after         deducing the respective capacitive currents from the original         currents.

Writing (21) in a more compact form results in:

$\begin{matrix} {{{{{\underset{\_}{A}}_{2}^{{comp\_}1}\left( d_{{comp\_}1} \right)}^{2} + {{\underset{\_}{A}}_{1}^{{comp\_}1}d_{{comp\_}1}} + {\underset{\_}{A}}_{0}^{{comp\_}1} + {{\underset{\_}{A}}_{00}^{{comp} - 1}R_{F}}} = 0}{{where}\text{:}}{{\underset{\_}{A}}_{2}^{{comp\_}1} = {{A_{2{\_ Re}}^{{comp\_}1} + {j\; A_{2{\_ Im}}^{{comp\_}1}}} = {{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{1L}^{long}}}}{{\underset{\_}{A}}_{1}^{{comp\_}1} = {{A_{1{\_ Re}}^{{comp\_}1} + {j\; A_{1{\_ Im}}^{{comp\_}1}}} = {{{\underset{\_}{L}}_{1}{\underset{\_}{Z}}_{1L}^{long}} - {{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{A\_ p}^{{comp\_}1}}}}}{{\underset{\_}{A}}_{0}^{{comp\_}1} = {{A_{0{\_ Re}}^{{comp\_}1} + {j\; A_{0{\_ Im}}^{{comp\_}1}}} = {{- {\underset{\_}{L}}_{1}}{\underset{\_}{Z}}_{A\_ p}^{{comp\_}1}}}}{{\underset{\_}{A}}_{00}^{{comp\_}1} = {{A_{00{\_ Re}}^{{comp\_}1} + {j\; A_{00{\_ Im}}^{{comp\_}1}}} = \frac{{\underset{\_}{M}}_{1}\left( {{{\underset{\_}{a}}_{F1}\Delta\;{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{AA2}}} \right)}{{\underset{\_}{I}}_{A\_ p}^{{comp\_}1}}}}} & \left( {21a} \right) \end{matrix}$

${\underset{\_}{Z}}_{A\_ p}^{{comp\_}1} = \frac{{\underset{\_}{V}}_{A\_ p}}{{\underset{\_}{I}}_{A\_ p}^{{comp\_}1}}$

-   -   fault loop impedance calculated with:         V_(A) _(—) _(p)—original (uncompensated) fault loop voltage (2),         I_(A) _(—) _(p) ^(comp) ^(—) ¹=a₁I_(A1) _(—) _(comp) _(—)         ₁+a₂I_(A2) _(—) _(comp) _(—) ₁+a₀I_(A0) _(—) _(comp) _(—)         ₁—fault loop current composed of the positive (12),         negative (16) and zero (17) sequence currents obtained after         deducing the respective capacitive currents from the original         currents,

K₁, L₁, M₁—coefficients gathered in TABLE 3.

Formula (21a) can be written down separately for real and imaginary parts: A ₂ _(—) _(Re) ^(comp) ^(—) ¹(d _(comp) _(—) ₁)² +A ₁ _(—) _(Re) ^(comp) ^(—) ¹ d _(comp) _(—) ₁ +A ₀ _(—) _(Re) ^(comp) ^(—) ¹ +A ₀₀ _(—) _(Re) ^(comp) ^(—) ¹ R _(F)=0  (21b) A ₂ _(—) _(Im) ^(comp) ^(—) ¹(d _(comp) _(—) ₁)² +A ₁ _(—) _(Im) ^(comp) ^(—) ¹ d _(comp) _(—) ₁ +A ₀ _(—) _(Im) ^(comp) ^(—) ¹ +A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ R _(F)=0  (21c)

Combining (21b) and (21c) in such the way that fault resistance is eliminated, that is, equation (21b) is multiplied by A₀₀ _(—) _(Im) ^(comp) ^(—) ¹ and equation (21c) by A₀₀ _(—) _(Re) ^(comp) ^(—) ¹ and then subtracting them, yields the quadratic formula for a sought fault distance: B ₂ ^(comp) ^(—) ¹(d _(comp) _(—) ₁)² +B ₁ ^(comp) ^(—) ¹ d _(comp) _(—) ₁ +B ₀ ^(comp) ^(—) ¹=0  (22) where: B ₂ ^(comp) ^(—) ¹ =A ₂ _(—) _(Re) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ −A ₂ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Re) ^(comp) ^(—) ¹ B ₁ ^(comp) ^(—) ¹ =A ₁ _(—) _(Re) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ −A ₁ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Re) ^(comp) ^(—) ¹ B ₀ ^(comp) ^(—) ¹ =A ₀ _(—) _(Re) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ −A ₀ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Re) ^(comp) ^(—) ¹

Equation (22) has two roots [(d_(comp) _(—) ₁)₁, (d_(comp) _(—) ₁)₂] for the distance to fault:

$\begin{matrix} {{\left( d_{{comp\_}1} \right)_{1} = \frac{{- B_{1}^{{comp\_}1}} - \sqrt{\left( B_{1}^{{comp\_}1} \right)^{2} - {4B_{2}^{{comp\_}1}B_{0}^{{comp\_}1}}}}{2B_{2}^{{comp\_}1}}}{\left( d_{{comp\_}1} \right)_{2} = \frac{{- B_{1}^{{comp\_}1}} + \sqrt{\left( B_{1}^{{comp\_}1} \right)^{2} - {4B_{2}^{{comp\_}1}B_{0}^{{comp\_}1}}}}{2B_{2}^{{comp\_}1}}}} & (23) \end{matrix}$

The root, which corresponds to the selected previously root (10) for d (uncompensated) is taken as the valid result. The compensation procedure requires iterative calculations, performed until the convergence is achieved (i.e. until the location estimates cease to change from the previous estimates) or as with a fixed number of iterations such as 2-3 iterations. The calculated distance to a fault from a particular (say, present iteration) is utilized for determining the shunt current in the next iteration.

The method of the invention is illustrated in two flow-charts of the FL algorithm, FIG. 10, single line and in FIG. 11, parallel lines.

As shown in the flowchart in FIG. 10 the following measurements are utilized:

-   -   voltages from the side A from particular phases a, b, c: v_(A)         _(—) _(a), v_(A) _(—) _(b), v_(A) _(—) _(c)     -   currents from the side A from particular phases a, b, c: i_(A)         _(—) _(a), i_(A) _(—) _(b), i_(A) _(—) _(c)

The input data utilized at step 101, Input data and measurements, are as follows:

-   -   impedances of the line for the positive (Z_(1L)) and zero         (Z_(0L)) sequences,     -   impedance of the extra link 25, 45, 55 between the substations         A, B for the positive (negative) sequences (Z_(1AB))     -   source impedances for the positive (negative) sequences         (Z_(1sA), Z_(1sB)): the representative values or the measured         values are used, and a communication means is used for sending         the measured remote source impedance as previously described,     -   information on the fault type (from the protective relay).

The measured fault quantities (voltages and currents) undergo adaptive filtering at step 104, Adaptive filtering of phase quantities, aimed at rejecting the dc components from currents and the transients induced by Capacitive Voltage Transformers (CVTs).

In the next step the symmetrical components of voltages and currents are calculated, step 105, which is equivalent to that as shown in FIGS. 3 a, 3 b. The fault loop signals are composed: the fault loop voltage as in (2), while the fault loop current as in (3)—but with taking a_(0m)=0.

The distance to fault without taking into account the shunt capacitances effect (d) is calculated at step 106 by solving the quadratic formula (9). The solution of (9) is presented in (10).

The result obtained without taking into account the shunt capacitances effect d following 106 is treated as the starting value for performing the compensation for shunt capacitances. The distributed long line model is applied for the compensation.

The following additional data is required for the calculating the compensation for shunt capacitance, step 107:

-   -   positive sequence capacitance of the line (C_(1L))     -   zero sequence capacitance of the line (C_(0L))     -   line length (1), which is used to express line         impedances/capacitances per km length.

The first iteration of the compensation leads to the quadratic equation (22), which is solved in (23). The next iterations are performed analogously. Iterative calculations are performed until the convergence is achieved or a fixed number of iterations, i.e. 2-3 iterations, may be made. The calculated distance to a fault from a particular (say, present iteration) is utilized for determining the shunt current in the next iteration. After completing the iterative calculations the distance to fault d_(comp) is obtained.

As shown in the flowchart in FIG. 11 for parallel lines the following measurements are utilized:

-   -   voltages from the side A and line LA from particular phases a,         b, c: v_(AA) _(—) _(a), v_(AA) _(—) _(b), v_(AA) _(—) _(c)     -   currents from the side A and line LA from particular phases a,         b, c: i_(AA) _(—) _(a), i_(AA) _(—) _(b), i_(AA) _(—) _(c)     -   zero sequence current from the healthy parallel line LB: i_(AB0)

The input data utilized in step 111, Input data and measurements, are as follows:

-   -   impedances of the faulted line for the positive (Z_(1LA)) and         zero (Z_(0LA)) sequences     -   impedance of the healthy line for the positive (negative)         sequence (Z_(1LB))     -   impedance of the extra link between the substations A, B for the         positive (negative) sequences (Z_(1AB))     -   zero sequence impedance for the mutual coupling (Z_(0m))     -   representative values of the source impedances for the positive         (negative) sequences (Z_(1sA), Z_(1sB))     -   information on the fault type is obtained from the protective         relay.

The measured fault quantities, the voltages and currents, undergo adaptive filtering at step 114 for the purpose of rejecting the dc components from currents and the transients induced by Capacitive Voltage Transformers (CVTs).

In the next step 115 the symmetrical components of voltages and currents are calculated as shown in FIGS. 3 a, 3 b. The fault loop signals are composed: the fault loop voltage as in (2), while the fault loop current as in (3)—but with taking a_(0m)=a₀.

The distance to fault without taking into account the shunt capacitances effect (d) is calculated in step 116 by solving the quadratic formula (9). The solution of (9) is presented in (10)

The result obtained without taking into account the shunt capacitances effect (d) is treated as the starting value for performing the compensation for shunt capacitances. The distributed long line model is applied for the compensation.

The following additional data for the faulted line is required for the compensation for shunt capacitances step 117:

-   -   positive sequence capacitance of the line (C_(1L))     -   zero sequence capacitance of the line (C_(0L))     -   line length (1), which is used to express line         impedances/capacitances per km length.

In the case of a single line, the compensation is performed analogously. Iterative calculations are performed until the convergence is achieved or by using a fixed number of iterations, e.g. 2-3 iterations. The calculated distance to a fault from a particular iteration, for example the present iteration, is utilized for determining the shunt current in the next iteration. After completing the iterative calculations the distance to fault d_(comp) is obtained.

FIG. 18 shows an embodiment of a device for determining the distance from one end, here shown as end A of a section of power transmission or distribution line A-B, to a fault F on the power line according to the method of the invention described. The fault locator device 1 receives measurements from measuring devices located at one end A such as current measuring means 14, and voltage measurements from voltage measurement means 11. The fault locator device may comprise measurement value converters, members for treatment of the calculating algorithms of the method, indicating means for the calculated distance to fault and a printer or connection to a printer or facsimile machine or similar for printout of the calculated fault. In a preferred embodiment of the device the fault locator comprises computer program means for providing a display of the information provided by the method of the invention, such as distance to a fault d or d_(comp) on a terminal which may be remote from the location of the line and/or the fault locator. Preferably the computer program means receives information from the fault locator and makes it available to provide information for a display of a computer such that an operator or engineer may see a value for the calculated distance to a fault displayed. The value may be displayed relative to a schematic display of the line or network in which the fault has taken place.

In the embodiment shown, measuring device 14 for continuous measurement of all the phase currents, and measuring device 11 for measurement of voltages are arranged at one end, station A. Optionally, measuring devices such as 15, 13 may also be arranged at station B but they are not necessary to practise the invention. The measured values such as: the three-phase voltages (V_(AA)) of the faulted line, three-phase voltages (I_(AA)) of the faulted line and zero sequence current (I_(AB0)) from the healthy parallel line (note that zero sequence is not present when the single line only is considered), and a value representative of the source impedance at B, Z_(1sB) as a are all passed to a calculating unit comprised in fault locator 1, filtered such as described in relation to FIGS. 3 a, 3 b, and stored in memory means. The calculating unit is provided with the calculating algorithms described, and programmed for the processes needed for calculating the distance to fault. Optionally, the source impedance for the remote end, Z_(1sB) may be measured by the remote device RD, 10, and the information sent via a high speed communication means 60 to the fault locator at A. In some applications it will be preferable to use a measured value sent from B instead of a representative value stored at A. It may be seen in FIG. 18 that current measuring means 15 and voltage measuring means 13 at remote end B may provide a RD 10, a fault locator, or any suitable device with measurements to calculate the remote source impedance

The calculating unit of fault locator 1 is provided with pre-fault phase currents and also with known values such as shunt capacitances and the impedances of the line. In respect of the occurrence of a fault, information regarding the type of fault, phase-to-phase, phase-to-ground etc., may be supplied to the calculating unit of the fault locator. When the calculating unit has determined the distance to fault, it is displayed on the device and/or sent to display means which may be remotely located. A printout or fax of the result may also be provided. In addition to signalling the fault distance, the device can produce reports, in which are recorded measured values of the currents of both lines, voltages, type of fault and other measured and/or calculated information associated with a given fault at a distance.

The method and a fault locator device according to any embodiment of the invention may be used to determine distance to a fault on a section of power transmission line. The present invention may also be used to determine a distance to a fault on a section of a power distribution line, or any other line or bus arranged for any of generation, transmission, distribution, control or consumption of electrical power.

The fault locator device and system may comprise filters for filtering the signals, converters for sampling the signals and one or more micro computers. The microprocessor (or processors) comprises a central processing unit CPU performing the steps of the method according to the invention. This is performed with the aid of a computer program, which is stored in the program memory. It is to be understood that the computer program may also be run on one or more general purpose industrial computers or microprocessors instead of a specially adapted computer.

The computer program comprises computer program code elements or software code portions that make the computer perform the method using equations, algorithms, data and calculations previously described. A part of the program may be stored in a processor as above, but also in a ROM, RAM, PROM or EPROM chip or similar. The program in part or in whole may also be stored on, or in, other suitable computer readable medium such as a magnetic disk, CD-ROM or DVD disk, hard disk, magneto-optical memory storage means, in volatile memory, in flash memory, as firmware, or stored on a data server.

A computer program according to the invention may be stored at least in part on different mediums that are computer readable. Archive copies may be stored on standard magnetic disks, hard drives, CD or DVD disks, or magnetic tape. The databases and libraries are stored preferably on one or more local or remote data servers, but the computer program and/or computer program product may, for example at different times, be stored in any of; a volatile Random Access memory (RAM) of a computer or processor, a hard drive, an optical or magneto-optical drive, or in a type of non-volatile memory such as a ROM, PROM, or EPROM device. The computer program may also be arranged in part as a distributed application capable of running on several different computers or computer systems at more or less the same time.

In another preferred embodiment, the fault locator may be used with parallel lines to locate single phase-to-ground faults (a-g, b-g, c-g faults) in case of unavailability of measurements from the healthy parallel line. Two modes of the healthy line operation are taken into account:

-   -   the healthy line being in operation,     -   the healthy line switched-off and grounded at both the ends.

FIG. 19 shows the flow-chart of the algorithm for locating faults in parallel lines under unavailability of measurements from the healthy parallel line. The unavailable, here healthy line zero sequence current, is required for reflecting the mutual coupling effect under single phase-to-ground faults (a-g, b-g, c-g faults). The unavailable current is thus estimated. The other faults can be located with the standard fault location algorithm (such as the algorithm from reference [1]).

The sequence of computations for the presented one-end fault location algorithm is as follows.

As shown in the flowchart in FIG. 19 the following measurements are utilized:

-   -   voltages from the side A and line LA (faulted) from particular         phases: v_(AA) _(—) _(a), v_(AA) _(—) _(b), v_(AA) _(—) _(c),     -   currents from the side A and line LA (faulted) from particular         phases: i_(AA) _(—) _(a), i_(AA) _(—) _(b), i_(AA) _(—) _(c)

The utilized input data are as follows:

-   -   impedances of the faulted line for the positive (Z_(1LA)) and         zero (Z_(0LA)) sequences     -   impedance of the healthy line for the zero sequence (Z_(0LB))     -   zero sequence impedance for the mutual coupling (Z_(0m))     -   information on the fault type (from the protective relay)     -   mode of the healthy line operation: in operation or switched-off         and grounded at both ends).

The measured fault quantities (voltages and currents) undergo adaptive filtering aimed at rejecting the dc components from currents and the transients induced by Capacitive Voltage Transformers (CVTs).

The following equations are used in the method of the present parallel line embodiment. In addition to the algorithm (1) described above for estimating distance to fault (d [pu]) by considering the Kirchhoff's voltage law for the fault loop as seen from the locator installation point:

$\begin{matrix} {{{\underset{\_}{V}}_{AA\_ p} - {d{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ p}} - {R_{F}\left( {\frac{{\underset{\_}{a}}_{F1}\Delta\;{\underset{\_}{I}}_{AA1}}{{\underset{\_}{k}}_{F1}} + \frac{{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{AA2}}{{\underset{\_}{k}}_{F2}} + \frac{{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{AA0}}{{\underset{\_}{k}}_{F0}}} \right)}} = 0} & (24) \end{matrix}$ the fault loop voltage (V_(AA) _(—) _(p)) and current (I_(AA) _(—) _(p)) can be expressed in terms of the symmetrical quantities:

$\begin{matrix} {{{\underset{\_}{V}}_{AA\_ p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{AA2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{AA0}}}}{and}} & (2) \\ {{{\underset{\_}{I}}_{AA\_ p} = {{\underset{\_}{I}}_{AA\_ p}^{SL} + {\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}I_{AB0}}}}{{where}\text{:}}} & (25) \\ {{\underset{\_}{I}}_{AA\_ p}^{SL} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{AA2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{AA0}}}} & \left( {25a} \right) \end{matrix}$ is the fault loop current without compensating for the mutual coupling effect (i.e. composed as for the single line—superscript SL), a₁, a₂, a₀—complex coefficients gathered in TABLE 1 (the derivation as in APPENDIX A1), V_(AA1), V_(AA2), V_(AA0)—positive, negative and zero sequence of measured voltages, I_(AA1), I_(AA2), I_(AA0)—positive, negative and zero sequence current from the faulted line LA, I_(AB0)—unavailable zero sequence current from the healthy parallel line LB (to be estimated), Z_(1LA), Z_(0LA)—positive and zero sequence impedance the whole line LA, Z_(0m)—zero sequence impedance for mutual coupling between the lines LA and LB, R_(F)—unknown fault resistance.

In the next step the symmetrical components of voltages and currents are calculated as shown in FIGS. 3 a, 3 b. The fault loop signals are composed: the fault loop voltage as in (2), while the fault loop current as in (25). The formulae (2)-(25) present the fault loop signals expressed in terms of the symmetrical components of the measured signals. However, one can use the classic way for composing the fault loop signals.

The presented method covers single phase-to ground faults (a-g, b-g, c-g faults). The other remaining faults have to be located with the fault location algorithms described above or standard fault location algorithm [1]. Distance to fault (d) for the considered here single phase-to ground faults is calculated by solving the quadratic formula for a sought distance to fault (26). Equation (26) is the same as equation (10) except that the values for B₁,B₂,B₃ are different to the values determined in (10). The solution gives two roots:

$\begin{matrix} {{d_{1} = \frac{{- B_{1}} - \sqrt{B_{1}^{2} - {4B_{2}B_{0}}}}{2B_{2}}}{d_{2} = \frac{{- B_{1}} + \sqrt{B_{1}^{2} - {4B_{2}B_{0}}}}{2B_{2}}}} & (26) \end{matrix}$ (as above, the root which fulfils the condition (0≦d≦1) is selected as the solution for the distance to fault). One has to substitute the following into (26): B ₂=real(A ₂)imag(A ₀₀)−real(A ₀₀)imag(A ₂) B ₁=real(A ₁)imag(A ₀₀)−real(A ₀₀)imag(A ₁) B ₀=real(A ₀)imag(A ₀₀)−real(A ₀₀)imag(A ₀) where:

${\underset{\_}{A}}_{2} = {{{- {\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{K}}_{1}{\underset{\_}{I}}_{AA\_ p}^{SL}} - {\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{P}}_{0}}{\underset{\_}{K}}_{1}{\underset{\_}{I}}_{{AA}\; 0}} - {\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{P}}_{0}}{\underset{\_}{Q}}_{0}}}$ ${\underset{\_}{A}}_{1} = {{{\underset{\_}{K}}_{1}{\underset{\_}{V}}_{AA\_ p}} - {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{L}}_{1}{\underset{\_}{I}}_{AA\_ p}^{SL}} - {\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{P}}_{0}}{\underset{\_}{L}}_{1}{\underset{\_}{I}}_{{AA}\; 0}}}$ ${{\underset{\_}{A}}_{0} = {{\underset{\_}{L}}_{1}{\underset{\_}{V}}_{AA\_ p}}}\;$ ${\underset{\_}{A}}_{00} = {{- {\underset{\_}{M}}_{1}}{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{{AA2}\;}}$ ${\underset{\_}{Q}}_{0} = {{\underset{\_}{M}}_{1}\left( {{{\underset{\_}{b}}_{F\; 1}\Delta\;{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{b}}_{F\; 2}{\underset{\_}{I}}_{{AA}\; 2}}} \right)}$

The recommended SET of the coefficients b_(F1), b_(F2) are taken from TABLE 4 and the recommended SET of the coefficients a_(F1), a_(F2), a_(F0) from TABLE 5.

Fault loop voltage in this embodiment is found from the TABLE below

Fault loop voltage composed Fault loop voltage composed in terms of symmetrical components as in the classic approach ${\underset{\_}{V}}_{{AA}_{-}p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{{AA}\; 2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{{AA}\; 0}}}$ a-g fault: a ₁ = a ₂ = a ₀ = 1 ${a\text{-}g\mspace{14mu}{fault}\text{:}\mspace{20mu}{\underset{\_}{V}}_{{AA}_{-}p}} = {\underset{\_}{V}}_{{AA}_{-}a}$ b-g fault: a ₁ = a ², a ₂ = a, a ₀ = 1 ${a\text{-}g\mspace{14mu}{fault}\text{:}\mspace{20mu}{\underset{\_}{V}}_{{AA}_{-}p}} = {\underset{\_}{V}}_{{AA}_{-}b}$ c-g fault: a ₁ = a, a ₂ = a ², a ₀ = 1 ${a\text{-}g\mspace{14mu}{fault}\text{:}\mspace{20mu}{\underset{\_}{V}}_{{AA}_{-}p}} = {\underset{\_}{V}}_{{AA}_{-}c}$ ${\underset{\_}{a} = {\exp\left( {j\; 2{\pi/3}} \right)}},\mspace{25mu}{j = \sqrt{- 1}}$

Fault loop current I_(AA) _(—) _(p) ^(SL) composed as for the single line is found from the TABLE below

$\begin{matrix} {{\underset{\_}{I}}_{AA\_ p}^{SL}\mspace{11mu}{in}\mspace{14mu}{terms}\mspace{14mu}{of}\mspace{14mu}{symmetrical}} \\ {components} \end{matrix}\quad$ ${\underset{\_}{I}}_{AA\_ p}^{SL}\mspace{11mu}{as}\mspace{14mu}{in}\mspace{14mu}{the}\mspace{14mu}{classic}\mspace{14mu}{approach}$ $\begin{matrix} {{\underset{\_}{I}}_{AA\_ p}^{SL} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{AA2}} +}} \\ {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}\mspace{11mu}{\underset{\_}{I}}_{AA0}} \end{matrix}\quad$ a-g fault: a ₁ = a ₂ = a ₀ = 1 ${a\text{-}g\mspace{14mu}{fault}\text{:}\mspace{20mu}{\underset{\_}{I}}_{AA\_ p}^{SL}} = {{\underset{\_}{I}}_{AA\_ a} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}}}$ b-g fault: a ₁ = a ², a ₂ = a, a ₀ = 1 ${a\text{-}g\mspace{14mu}{fault}\text{:}\mspace{20mu}{\underset{\_}{I}}_{AA\_ p}^{SL}} = {{\underset{\_}{I}}_{AA\_ a} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}}}$ c-g fault: a ₁ = a, a ₂ = a ², a ₀ = 1 ${a\text{-}g\mspace{14mu}{fault}\text{:}\mspace{20mu}{\underset{\_}{I}}_{AA\_ p}^{SL}} = {{\underset{\_}{I}}_{AA\_ a} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}}}$ ${\underset{\_}{a} = {\exp\left( {j\; 2{\pi/3}} \right)}},\mspace{25mu}{j = \sqrt{- 1}}$ ${{where}\text{:}\mspace{20mu}{\underset{\_}{k}}_{0}} = \frac{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{1{LA}}}{{\underset{\_}{Z}}_{1{LA}}}$

The complex coefficients dependent on the mode of the healthy parallel line operation:

a) healthy line LB is in operation:

${\underset{\_}{P}}_{0} = \frac{{\underset{\_}{Z}}_{0{LB}} - {\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{0m}}$ ${\underset{\_}{K}}_{1} = {- {{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}} + {\underset{\_}{Z}}_{1{LB}}} \right)}}$ ${\underset{\_}{L}}_{1} = {{- {\underset{\_}{K}}_{1}} + {{\underset{\_}{Z}}_{1{LB}}{\underset{\_}{Z}}_{1{sB}}}}$ ${\underset{\_}{M}}_{1} = {{{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{Z}}_{1{LB}}} + {{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}}} \right)} + {{\underset{\_}{Z}}_{1{LB}}\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}}} \right)}}$ b) healthy line LB is switched-off and grounded:

${\underset{\_}{P}}_{0} = {- \frac{{\underset{\_}{Z}}_{0{LB}}}{{\underset{\_}{Z}}_{0m}}}$ ${\underset{\_}{K}}_{1} = {- {\underset{\_}{Z}}_{1{LA}}}$ ${\underset{\_}{L}}_{1} = {{\underset{\_}{Z}}_{1{LA}} + {\underset{\_}{Z}}_{1{sB}}}$ ${\underset{\_}{M}}_{1} = {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{LA}}}$

In another embodiment, which is presented here, the method is applied for the standard availability of the measured signals and is only valid for ground faults, including both:

-   -   single phase-to-ground faults     -   phase-to-phase-to-ground faults.

Fault location procedures obtained for these faults under the standard availability of the fault locator input signals are extremely simple and compact. Distance to a fault is calculated with a first order formula.

FIG. 22 shows a flow-chart of the developed algorithm for locating ground faults in parallel transmission lines. The sequence of computations for the presented one-end fault locator is as follows. As shown in the flowchart in FIG. 22.B3 the following measurements are utilized:

-   -   voltages from the side A and line LA from particular phases a,         b, c: v_(AA) _(—) _(a), v_(AA) _(—) _(b), v_(AA) _(—) _(c)     -   currents from the side A and line LA from particular phases a,         b, c: i_(AA) _(—) _(a), i_(AA) _(—) _(b), i_(AA) _(—) _(c)     -   zero sequence current from the healthy parallel line LB: I_(AB0)

The utilized input data are as follows:

-   -   impedances of the faulted line for the positive (Z_(1LA)) and         zero (Z_(0LA)) sequences,     -   impedance of the healthy line for the zero sequence (Z_(0LB))     -   zero sequence impedance for the mutual coupling (Z_(0m))     -   information on the fault type from the protective relay.

The measured fault quantities (voltages and currents) undergo adaptive filtering aimed at rejecting the dc components from currents and the transients induced by Capacitive Voltage Transformers (CVTs), preferably as described and shown in relation to FIGS. 3 a, 3 b.

A generalized model of the fault loop, used for deriving the algorithm of the present embodiment, is stated as follows: V _(AA) _(—) _(p) −dZ _(1LA) I _(AA) _(—) _(p) −R _(F)(a _(F1) I _(F1) +a _(F2) I _(F2) +a _(F0) I _(F0))=0  (1) where: d—unknown and sought distance to fault, Z_(1LA)—positive sequence impedance of the faulted line, V_(AA) _(—) _(p), I_(AA) _(—) _(p)—fault loop voltage and current composed according to the fault type, R_(F)—fault resistance, I_(Fi)—sequence components of the total fault current (i=0—zero sequence, i=1 positive sequence, i=2—negative sequence), a_(Fi)—weighting coefficients (TABLE 2).

Fault loop voltage and current can be expressed as in classic distance protection technique or, as in this document, in terms of the local measurements and with using the coefficients (a₀, a₁, a₂) which are gathered in TABLE 1 (derivation of the coefficients is contained in APPENDIX APP.1): V _(AA) _(—) _(p) =a ₁ V _(AA1) +a ₂ V _(AA2) +a ₀ V _(AA0)  (2)

$\begin{matrix} {{\underset{\_}{I}}_{AA\_ p} = \;{{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{{AA}\; 2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AA}\; 0}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AB}\; 0}}}} & (3) \end{matrix}$ where: AA, AB—subscripts used for indicating measurements acquired from the faulted line (AA) and from the healthy line (AB), respectively, Z_(0LA), Z_(0m)—impedance of the faulted line and mutual coupling between the lines for the zero sequence, respectively.

In the next step the symmetrical components of voltages and currents are calculated as shown in FIGS. 3 a, 3 b. The fault loop signals are composed: the fault loop voltage as in (2), while the fault loop current as in (3). The formulae (2)-(3) present the fault loop signals expressed in terms of the symmetrical components of the measured signals. However, one may instead use the classic way for composing the fault loop signals, as shown in APPENDIX A1.

The presented method covers single phase-to ground faults (a-g, b-g, c-g faults) and phase-to-phase-to-ground faults (a-b-g, b-c-g, c-a-g faults), thus, the faults for which the highest fault resistance can be expected. The other remaining faults have to be located with the fault location algorithms described above or a standard fault location algorithm, such as for example the fault locator from reference [1]).

Distance to fault (d) for single phase-to-ground faults is calculated as follows:

a-g fault:

$\begin{matrix} {d = \frac{{imag}\left\{ {{\underset{\_}{V}}_{AA\_ p}\left\lbrack {3\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}{{imag}\left\{ {\left( {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ P}} \right)\left\lbrack {3\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}} & \left( {27a} \right) \end{matrix}$ where:

Fault loop signals composed Fault loop signals composed in terms of symmetrical components as in the classic approach ${\underset{\_}{V}}_{{AA}_{-}p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{{AA}\; 2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{{AA}\; 0}}}$ ${\underset{\_}{V}}_{{AA}_{-}p} = {\underset{\_}{V}}_{{AA}_{-}a}$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{{AA}\; 2}} +}} \\ {{{+ {\underset{\_}{a}}_{0}}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AA}\; 0}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AB}\; 0}}} \end{matrix}\quad$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{\underset{\_}{I}}_{AA\_ a} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}} +}} \\ {{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{{AB}\; 0}} \end{matrix}\quad$ where: a ₁ = a ₂ = a ₀ = 1 $\begin{matrix} {{{{where}\text{:}\mspace{20mu}{\underset{\_}{k}}_{0}} = \frac{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{1{LA}}}{{\underset{\_}{Z}}_{1{LA}}}},} \\ {{\underset{\_}{k}}_{0m} = \frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1\;{LA}}}} \end{matrix}\quad$

${\underset{\_}{P}}_{0} = {{\frac{{\underset{\_}{Z}}_{0{LB}} - {\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{0m}}\mspace{14mu}\text{(for~~symmetrical~~lines:}\mspace{14mu}{\underset{\_}{P}}_{0}} = {1{\text{)}.}}}$ b-g fault:

$\begin{matrix} {d = \frac{{imag}\left\{ {{\underset{\_}{V}}_{AA\_ p}\left\lbrack {3{{\underset{\_}{a}}^{2}\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)}} \right\rbrack}^{*} \right\}}{{imag}\left\{ {\left( {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ P}} \right)\left\lbrack {3{{\underset{\_}{a}}^{2}\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)}} \right\rbrack}^{*} \right\}}} & \left( {27b} \right) \end{matrix}$ where:

Fault loop signals composed Fault loop signals composed in terms of symmetrical components as in the classic approach ${\underset{\_}{V}}_{{AA}_{-}p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{{AA}\; 2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{{AA}\; 0}}}$ ${\underset{\_}{V}}_{{AA}_{-}p} = {\underset{\_}{V}}_{{AA}_{-}b}$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{{AA}\; 2}} +}} \\ {{{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AA}\; 0}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AB}\; 0}}} \end{matrix}\quad$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{\underset{\_}{I}}_{AA\_ b} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}} +}} \\ {{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{{AB}\; 0}} \end{matrix}\quad$ where: a ₁ = a ², a ₂ = a, a ₀ = 1a = exp(j2π/3) $\begin{matrix} {{{{where}\text{:}\mspace{20mu}{\underset{\_}{k}}_{0}} = \frac{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{1{LA}}}{{\underset{\_}{Z}}_{1{LA}}}},} \\ {{\underset{\_}{k}}_{0m} = \frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1\;{LA}}}} \end{matrix}\quad$

${\underset{\_}{P}}_{0} = {{\frac{{\underset{\_}{Z}}_{0{LB}} - {\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{0m}}\mspace{14mu}\text{(for~~symmetrical~~lines:}\mspace{14mu}{\underset{\_}{P}}_{0}} = {1{\text{)}.}}}$ c-g fault:

$\begin{matrix} {d = \frac{{imag}\left\{ {{\underset{\_}{V}}_{AA\_ p}\left\lbrack {3{\underset{\_}{a}\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)}} \right\rbrack}^{*} \right\}}{{imag}\left\{ {\left( {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ P}} \right)\left\lbrack {3{\underset{\_}{a}\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)}} \right\rbrack}^{*} \right\}}} & \left( {27c} \right) \end{matrix}$ where:

Fault loop signals composed Fault loop signals composed in terms of symmetrical components as in the classic approach ${\underset{\_}{V}}_{{AA}_{-}p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{{AA}\; 2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{{AA}\; 0}}}$ ${\underset{\_}{V}}_{{AA}_{-}p} = {\underset{\_}{V}}_{{AA}_{-}c}$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{{AA}\; 2}} +}} \\ {{{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AA}\; 0}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AB}\; 0}}} \end{matrix}\quad$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{\underset{\_}{I}}_{AA\_ c} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}} +}} \\ {{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{{AB}\; 0}} \end{matrix}\quad$ where: a ₁ = a, a ₂ = a ², a ₀ = 1a = exp(j2π/3) $\begin{matrix} {{{{where}\text{:}\mspace{20mu}{\underset{\_}{k}}_{0}} = \frac{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{1{LA}}}{{\underset{\_}{Z}}_{1{LA}}}},} \\ {{\underset{\_}{k}}_{0m} = \frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1\;{LA}}}} \end{matrix}\quad$

${\underset{\_}{P}}_{0} = {{\frac{{\underset{\_}{Z}}_{0{LB}} - {\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{0m}}\mspace{14mu}\text{(for~~symmetrical~~lines:}\mspace{14mu}{\underset{\_}{P}}_{0}} = {1{\text{)}.}}}$

Distance to fault (d) for phase-to- phase-to ground faults can be calculated in two different ways, depending whether the pre-fault currents can be used or have to be avoided.

1. Procedure for distance to fault calculation with the use of pre-fault measurements:

a-b-g fault:

$\begin{matrix} {d = \frac{{imag}\left\{ {{\underset{\_}{V}}_{AA\_ p}\left\lbrack {\underset{\_}{W}\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}{{imag}\left\{ {\left( {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ P}} \right)\left\lbrack {\underset{\_}{W}\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}} & \left( {28a} \right) \end{matrix}$

where: ${\underset{\_}{W} = \frac{{\underset{\_}{A}\Delta\;{\underset{\_}{I}}_{{AA}\; 1}} + {\underset{\_}{B\; I}}_{{AA}\; 2}}{{\underset{\_}{C}\Delta\;{\underset{\_}{I}}_{{AA}\; 1}} + {\underset{\_}{D\; I}}_{{AA}\; 2}}},{\underset{\_}{A} = {1 - {\underset{\_}{a}}^{2}}},{\underset{\_}{B} = {1 - \underset{\_}{a}}},{\underset{\_}{C} = {1 + {\underset{\_}{a}}^{2}}},{\underset{\_}{D} = {1 + \underset{\_}{a}}}$

Fault loop signals composed Fault loop signals composed in terms of symmetrical components as in the classic approach ${\underset{\_}{V}}_{{AA}_{-}p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{{AA}\; 2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{{AA}\; 0}}}$ ${\underset{\_}{V}}_{{AA}_{-}p} = {{\underset{\_}{V}}_{{AA}_{-}a} - {\underset{\_}{V}}_{{AA}_{-}b}}$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{{AA}\; 2}} +}} \\ {{{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AA}\; 0}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AB}\; 0}}} \end{matrix}\quad$ ${\underset{\_}{I}}_{{AA}_{-}p} = {{\underset{\_}{I}}_{{AA}_{-}a} - {\underset{\_}{I}}_{{AA}_{-}b}}$ where: a ₁ = 1 − a ², a ₂ = 1 − a, a ₀ = 0 a = exp(j2π/3)

${\underset{\_}{P}}_{0} = {{\frac{{\underset{\_}{Z}}_{0{LB}} - {\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{0m}}\mspace{14mu}\text{(for~~symmetrical~~lines:}\mspace{14mu}{\underset{\_}{P}}_{0}} = {1{\text{)}.}}}$ b-c-g fault:

$\begin{matrix} {{d = \frac{{imag}\left\{ {{\underset{\_}{V}}_{AA\_ p}\left\lbrack {\underset{\_}{W}\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}{{imag}\left\{ {\left( {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ P}} \right)\left\lbrack {\underset{\_}{W}\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}}{{where}\text{:}}{{\underset{\_}{W} = \frac{{\underset{\_}{A}\Delta\;{\underset{\_}{I}}_{{AA}\; 1}} + {\underset{\_}{B\; I}}_{{AA}\; 2}}{{\underset{\_}{C}\Delta\;{\underset{\_}{I}}_{{AA}\; 1}} + {\underset{\_}{D\; I}}_{{AA}\; 2}}},{\underset{\_}{A} = {{\underset{\_}{a}}^{2} - \underset{\_}{a}}},{\underset{\_}{B} = {\underset{\_}{a} - {\underset{\_}{a}}^{2}}},{\underset{\_}{C} = {- 1}},{\underset{\_}{D} = {- 1}}}} & \left( {28b} \right) \end{matrix}$

Fault loop signals composed Fault loop signals composed in terms of symmetrical components as in the classic approach ${\underset{\_}{V}}_{{AA}_{-}p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{{AA}\; 2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{{AA}\; 0}}}$ ${\underset{\_}{V}}_{{AA}_{-}p} = {{\underset{\_}{V}}_{{AA}_{-}b} - {\underset{\_}{V}}_{{AA}_{-}c}}$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{{AA}\; 2}} +}} \\ {{{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AA}\; 0}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AB}\; 0}}} \end{matrix}\quad$ ${\underset{\_}{I}}_{{AA}_{-}p} = {{\underset{\_}{I}}_{{AA}_{-}b} - {\underset{\_}{I}}_{{AA}_{-}c}}$ where: a ₁ = a ² − a, a ₂ = a − a ², a ₀ = 0 a = exp(j2π/3)

${\underset{\_}{P}}_{0} = {{\frac{{\underset{\_}{Z}}_{0{LB}} - {\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{0m}}\mspace{14mu}\text{(for~~symmetrical~~lines:}\mspace{14mu}{\underset{\_}{P}}_{0}} = {1{\text{)}.}}}$ c-a-g fault:

$\begin{matrix} {{d = \frac{{imag}\left\{ {{\underset{\_}{V}}_{AA\_ p}\left\lbrack {\underset{\_}{W}\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}{{imag}\left\{ {\left( {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ P}} \right)\left\lbrack {\underset{\_}{W}\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}}{{where}\text{:}}{{\underset{\_}{W} = \frac{{\underset{\_}{A}\Delta\;{\underset{\_}{I}}_{{AA}\; 1}} + {\underset{\_}{B\; I}}_{{AA}\; 2}}{{\underset{\_}{C}\Delta\;{\underset{\_}{I}}_{{AA}\; 1}} + {\underset{\_}{D\; I}}_{{AA}\; 2}}},{\underset{\_}{A} = {\underset{\_}{a} - 1}},{\underset{\_}{B} = {{\underset{\_}{a}}^{2} - 1}},{\underset{\_}{C} = {\underset{\_}{a} + 1}},{\underset{\_}{D} = {{\underset{\_}{a}}^{2} + 1}}}} & \left( {28c} \right) \end{matrix}$

Fault loop signals composed Fault loop signals composed in terms of symmetrical components as in the classic approach ${\underset{\_}{V}}_{{AA}_{-}p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{{AA}\; 2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{{AA}\; 0}}}$ ${\underset{\_}{V}}_{{AA}_{-}p} = {{\underset{\_}{V}}_{{AA}_{-}\; c} - {\underset{\_}{V}}_{{AA}_{-}a}}$ $\begin{matrix} {{\underset{\_}{I}}_{{AA\_}\; p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{{AA}\; 2}} +}} \\ {{{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AA}\; 0}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{{AB}\; 0}}} \end{matrix}\quad$ ${\underset{\_}{I}}_{{AA}_{-}p} = {{\underset{\_}{I}}_{{AA}_{-}c} - {\underset{\_}{I}}_{{AA}_{-}a}}$ where: a ₁ = a − 1, a ₂ = a ² − 1, a ₀ = 0 a = exp(j2π/3)

${\underset{\_}{P}}_{0} = {\frac{{\underset{\_}{Z}}_{0{LB}} - {\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{0m}}\mspace{14mu}{\left( {{{for}\mspace{14mu}{symmetrical}\mspace{14mu}{lines}\text{:}\mspace{14mu}{\underset{\_}{P}}_{0}} = 1} \right).}}$ (for symmetrical lines: P₀=1). 2. Procedure for distance to fault calculation without the use of pre-fault measurements: a-b-g fault:

$\begin{matrix} \begin{matrix} {d = \frac{{imag}\left\lbrack {\left( {{\underset{\_}{V}}_{a} + {\underset{\_}{V}}_{b}} \right)\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)^{*}} \right\rbrack}{{imag}\left\lbrack {{{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{I}}_{a} + {\underset{\_}{I}}_{b} + {2{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}} + {2{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{\;{{AB}\; 0}}}} \right)}\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)^{*}} \right\rbrack}} & \; \end{matrix} & \left( {29a} \right) \end{matrix}$ b-c-g fault:

$\begin{matrix} {d = \frac{{imag}\left\lbrack {\left( {{\underset{\_}{V}}_{b} + {\underset{\_}{V}}_{c}} \right)\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)^{*}} \right\rbrack}{{imag}\left\lbrack {{{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{I}}_{b} + {\underset{\_}{I}}_{c} + {2{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}} + {2{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{\;{{AB}\; 0}}}} \right)}\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)^{*}} \right\rbrack}} & \left( {29b} \right) \end{matrix}$ c-a-g fault:

$\begin{matrix} {{d = \frac{{imag}\left\lbrack {\left( {{\underset{\_}{V}}_{c} + {\underset{\_}{V}}_{a}} \right)\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)^{*}} \right\rbrack}{\begin{matrix} {{imag}\left\lbrack {{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{I}}_{c} + {\underset{\_}{I}}_{a} + {2{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}} +} \right.} \right.} \\ \left. {\left. {2{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{\;{{AB}\; 0}}} \right)\left( {{\underset{\_}{I}}_{AA0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)^{*}} \right\rbrack \end{matrix}}}\text{where:}{{{\underset{\_}{k}}_{0} = \frac{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{1{LA}}}{{\underset{\_}{Z}}_{1{LA}}}},{{\underset{\_}{k}}_{0m} = {\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}.}}}} & \left( {29c} \right) \end{matrix}$

It is also noted that while the above describes exemplifying embodiments of the invention, there are several variations and modifications which may be made to the disclosed solution without departing from the scope of the present invention as defined in the appended claims.

Tables

TABLE 1 Coefficients for determining the fault loop voltage (V _(AA) _(—) _(p) ) and current (I _(AA) _(—) _(p) ) in terms of symmetrical components as defined in (2) and (3). Fault type a ₁ a ₂ a ₀ a-g 1 1 1 b-g a ² a 1 c-g a a ² 1 a-b, a-b-g 1 − a ² 1 − a ² 0 a-b-c, a-b-c-g b-c, b-c-g a ² − a a − a ² 0 c-a, c-a-g a − 1 a ² − 1 0 ${\underset{\_}{a} = {\exp\left( {j\; 2{\pi/3}} \right)}},\mspace{25mu}{j = \sqrt{- 1}}$

TABLE 1A $\begin{matrix} {{Fault}\mspace{14mu}{loop}\mspace{14mu}{voltage}\mspace{14mu}\left( {\underset{\_}{V}}_{AA\_ p} \right)\mspace{20mu}{and}\mspace{14mu}{current}\mspace{14mu}\left( {\underset{\_}{I}}_{AA\_ p} \right)\mspace{14mu}{expressed}} \\ {{with}\mspace{14mu}{using}\mspace{14mu}{the}\mspace{14mu}{classic}\mspace{14mu}{{approach}.}} \end{matrix}\quad$ Fault type ${\underset{\_}{V}}_{AA\_ p}$ ${\underset{\_}{I}}_{AA\_ p}$ a-g ${\underset{\_}{V}}_{AA\_ a}$ ${\underset{\_}{I}}_{AA\_ a} + {{\underset{\_}{k}}_{0}I_{AA0}} + {{\underset{\_}{k}}_{0m}I_{AB0}}$ b-g ${\underset{\_}{V}}_{AA\_ b}$ ${\underset{\_}{I}}_{AA\_ b} + {{\underset{\_}{k}}_{0}I_{AA0}} + {{\underset{\_}{k}}_{0m}I_{AB0}}$ c-g ${\underset{\_}{V}}_{AA\_ c}$ ${\underset{\_}{I}}_{AA\_ c} + {{\underset{\_}{k}}_{0}I_{AA0}} + {{\underset{\_}{k}}_{0m}I_{AB0}}$ a-b, a-b-ga-b-c,a-b-c-g ${\underset{\_}{V}}_{AA\_ a} - {\underset{\_}{V}}_{AA\_ b}$ ${\underset{\_}{I}}_{AA\_ a} - {\underset{\_}{I}}_{AA\_ b}$ b-c, b-c-g ${\underset{\_}{V}}_{AA\_ b} - {\underset{\_}{V}}_{AA\_ c}$ ${\underset{\_}{I}}_{AA\_ b} - {\underset{\_}{I}}_{AA\_ c}$ c-a, c-a-g ${\underset{\_}{V}}_{AA\_ c} - {\underset{\_}{V}}_{AA\_ a}$ ${\underset{\_}{I}}_{AA\_ c} - {\underset{\_}{I}}_{AA\_ a}$ $\begin{matrix} {{{\underset{\_}{k}}_{0} = \frac{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{1{LA}}}{{\underset{\_}{Z}}_{1{LA}}}},{{\underset{\_}{k}}_{0m} = \frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}},} \\ {{\underset{\_}{a} = {\exp\left( {{j2\pi}/3} \right)}},{j = \sqrt{- 1}}} \end{matrix}\quad$

TABLE 2 Alternative sets of the weighting coefficients for determining a voltage drop across the fault path resistance Set I Fault (recommended) Set II Set III type a _(F1) a _(F2) a _(F0) a _(F1) a _(F2) a _(F0) a _(F1) a _(F2) a _(F0) a-g 0 3 0 3 0 0 1.5 1.5 0 b-g 0 3a 0 3a ² 0 0 1.5a ² 1.5a 0 c-g 0 3a ² 0 3a 0 0 1.5a 1.5a ² 0 a-b 0 1 − a 0 1 − a ² 0 0 0.5(1 − a ²) 0.5(1 − a) 0 b-c 0 a − a ² 0 a ² − a 0 0 0.5(a ² − a) 0.5(a − a ²) 0 c-a 0 a ² − 1 0 a − 1 0 0 0.5(a − 1) 0.5(a ² − 1) 0 a-b-g 1 − a ² 1 − a 0 1 − a ² 1 − a 0 1 − a ² 1 − a 0 b-c-g a ² − a a − a ² 0 a ² − a a − a ² 0 a ² − a a − a ² 0 c-a-g a − 1 a ² − 1 0 a − 1 a ² − 1 0 a − 1 a ² − 1 0 a-b-c-g 1 − a ² 0 0 1 − a ² 0 0 1 − a ² 0 0 (a-b-c) a = exp(j2π/3), j = {square root over (−1)}

TABLE 3 Coefficients for determining a fault current distribution factor (6) SINGLE LINE (FIG. 4) ${\underset{\_}{Z}}_{1{AB}} \neq \propto$ $\begin{matrix} {{\underset{\_}{K}}_{1} = {{{- {\underset{\_}{Z}}_{1L}}{\underset{\_}{Z}}_{1{AB}}} - {\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}}} \right){\underset{\_}{Z}}_{1L}}}} \\ {{\underset{\_}{L}}_{1} = {{{\underset{\_}{Z}}_{1L}\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}}} \right)} + {{\underset{\_}{Z}}_{1{AB}}\left( {{\underset{\_}{Z}}_{1L} + {\underset{\_}{Z}}_{1{sB}}} \right)}}} \\ {{\underset{\_}{M}}_{1} = {{\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}}} \right)\left( {{\underset{\_}{Z}}_{1{AB}} + {\underset{\_}{Z}}_{1L}} \right)} + {{\underset{\_}{Z}}_{1L}{\underset{\_}{Z}}_{1{AB}}}}} \end{matrix}\quad$ $\left. {\underset{\_}{Z}}_{1{AB}}\rightarrow \propto \right.$ $\begin{matrix} {{\underset{\_}{K}}_{1} = {- {\underset{\_}{Z}}_{1L}}} \\ {{\underset{\_}{L}}_{1} = {{\underset{\_}{Z}}_{1L} + {\underset{\_}{Z}}_{1{sB}}}} \\ {{\underset{\_}{M}}_{1} = {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}} + {\underset{\_}{Z}}_{1L}}} \end{matrix}\quad$ PARALLEL LINES (FIG. 5) $\begin{matrix} {{\underset{\_}{K}}_{1} = {- {{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}} + {\underset{\_}{Z}}_{{1{LB}} + {AB}}} \right)}}} \\ {{\underset{\_}{L}}_{1} = {{{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}} + {\underset{\_}{Z}}_{{{1{LB}}\&}{AB}}} \right)} + {{\underset{\_}{Z}}_{{{1{LB}}\&}{AB}}{\underset{\_}{Z}}_{1{sB1}}}}} \\ \begin{matrix} {{\underset{\_}{M}}_{1} = {{{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{Z}}_{{{1{LB}}\&}{AB}}} + {{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}}} \right)} +}} \\ {{\underset{\_}{Z}}_{{{1{LB}}\&}{AB}}\left( {{\underset{\_}{Z}}_{1{sA}} + {\underset{\_}{Z}}_{1{sB}}} \right)} \end{matrix} \end{matrix}\quad$

TABLE 4 The recommended SET of the coefficients b_(F1), b_(F2) in relation to (26) Fault b_(F1) b_(F2) a-g 0 1 b-g 0 a² c-g 0 a a = exp(j2π/3), j = {square root over (−1)}

TABLE 5 The recommended SET of the coefficients a_(F1), a_(F2), a_(F0) in relation to (26) FAULT a_(F1) a_(F2) a_(F0) a-g 0 3 0 b-g 0 3a 0 c-g 0 3a² 0 a = exp(j2π/3), j = {square root over (−1)}

REFERENCES

-   [1] Eriksson L., Saha M. M., Rockefeller G. D.: An accurate fault     locator with compensation for apparent reactance in the fault     resistance resulting from remote-end infeed, IEEE Transactions on     PAS, Vol. PAS-104, No. 2, February 1985, pp. 424-436. -   [2] Saha M. M.: Method and device for locating a fault point on a     three-phase power transmission line. U.S. Pat. No. 4,559,491, Date     of Patent: Dec. 17, 1985. -   [3] McLAREN P. G., SWIFT G. W., ZHANG Z., DIRKS E., JAYASINGHLE R.     P., FERNANDO I., A new positive sequence directional element for     numerical distance relays, Proceedings of the Stockholm Power Tech     Conference, Stockholm, Sweden, 1995, pp. 540-545. -   [4] SAHA N. M., IZYKOWSKI J., KASZTENNY B., ROSOLOWSKI E., PALKI B.     S., Relaying algorithms for protection of series-compensated lines,     Proceedings of the International Conference on Modern Trends in the     Protection Schemes of Electric Power Apparatus and Systems, Oct.     28-30, 1998, New Delhi, India, pp. V-50-61. -   [5] NOVOSEL D., HART D. G., UDREN E., GARITTY J., Unsynchronized     two-terminal fault location estimation, IEEE Transactions on Power     Delivery, Vol. 11, No. 1, January 1996, pp. 130-138.

APPENDICES

A1. Derivation of the Coefficients a₁, a₂, a₀ (TABLE 1)

Single phase-to-ground fault: a-g fault V _(AA) _(—) _(p) =V _(AA) _(—) _(a) =V _(AA1) +V _(AA2) +V _(AA0) =a ₁ V _(AA1) +a ₂ V _(AA2) +a ₀ V _(AA0)

${\underset{\_}{I}}_{AA\_ p} = {{{\underset{\_}{I}}_{AA\_ a} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{AA0}} + {{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{AB0}}} = {{{\underset{\_}{I}}_{AA1} + {\underset{\_}{I}}_{AA2} + {\underset{\_}{I}}_{AA0} + {\frac{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{1{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{AA0}} + {\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{AB0}}} = {{{\underset{\_}{I}}_{AA1} + {\underset{\_}{I}}_{AA2} + {\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{AA0}} + {\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{AB0}}} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{AA2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0{LA}}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{AA0}} + {{\underset{\_}{a}}_{0m}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1{LA}}}{\underset{\_}{I}}_{AB0}}}}}}$

Thus: a₁=a₂=a₀=1

Inter-phase faults: a-b, a-b-g, a-b-c, a-b-c-g faults

${\underset{\_}{V}}_{AA\_ p} = {{{\underset{\_}{V}}_{AA\_ a} - {\underset{\_}{V}}_{AA\_ b}} = {{\left( {{\underset{\_}{V}}_{AA1} + {\underset{\_}{V}}_{AA2} + {\underset{\_}{V}}_{AA0}} \right) - \left( {{{\underset{\_}{a}}^{2}{\underset{\_}{V}}_{AA1}} + {\underset{\_}{aV}}_{AA2} + {\underset{\_}{V}}_{AA0}} \right)} = {{{\left( {1 - {\underset{\_}{a}}^{2}} \right){\underset{\_}{V}}_{AA1}} + {\left( {1 - \underset{\_}{a}} \right){\underset{\_}{V}}_{AA2}}} = {{{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{AA2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{AA0}{\underset{\_}{I}}_{AA\_ p}}} = {{{\underset{\_}{I}}_{AA\_ a} - {\underset{\_}{I}}_{AA\_ b}} = {{\left( {{\underset{\_}{I}}_{AA1} + {\underset{\_}{I}}_{AA2} + {\underset{\_}{I}}_{AA0}} \right) - \left( {{{\underset{\_}{a}}^{2}{\underset{\_}{I}}_{AA1}} + {\underset{\_}{aI}}_{AA2} + {\underset{\_}{I}}_{AA0}} \right)} = {{{\left( {1 - {\underset{\_}{a}}^{2}} \right){\underset{\_}{I}}_{AA1}} + {\left( {1 - \underset{\_}{a}} \right){\underset{\_}{I}}_{AA2}}} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{AA2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{L0}}{{\underset{\_}{Z}}_{L1}}{\underset{\_}{I}}_{AA0}}}}}}}}}}$

Thus: a₁=1-a², a₂=1-a, a₀=0

A2. Derivation of Coefficients a_(F1), a_(F2), a_(F0) (TABLE 2)

TABLE 2 contains three alternative sets (Set I, Set II, Set III) of the weighting coefficients, which are used for determining a voltage drop across a fault path. The coefficients are calculated from the boundary conditions—relevant for a particular fault type. This is distinctive that in all the sets the zero sequence is omitted (a_(F0)=0). This is advantage since the zero sequence impedance of a line is considered as the uncertain parameter. By setting a_(F0)=0 we limit adverse influence of the uncertainty with respect to the zero sequence impedance data upon the fault location accuracy. To be precise one has to note that this limitation is of course partial since it is related only to determining the voltage drop across a fault path. In contrast, while determining the voltage drop across a faulted line segment the zero sequence impedance of the line is used.

a-g fault, FIG. 12:

Taking into account that in the healthy phases: I_(F) _(—) _(b)=I_(F) _(—) _(c)=0 what gives:

${\underset{\_}{I}}_{F1} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {\underset{\_}{a\; I}}_{F\_ b} + {{\underset{\_}{a}}^{2}{\underset{\_}{I}}_{F\_ c}}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {\underset{\_}{a}0} + {{\underset{\_}{a}}^{2}0}} \right)} = {\frac{1}{3}{\underset{\_}{I}}_{F\_ a}}}}$ ${\underset{\_}{I}}_{F2} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {{\underset{\_}{a}}^{2}{\underset{\_}{I}}_{F\_ b}} + {\underset{\_}{a\; I}}_{F\_ c}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {{\underset{\_}{a}}^{2}0} + {\underset{\_}{a}0}} \right)} = {\frac{1}{3}{\underset{\_}{I}}_{F\_ a}}}}$ ${\underset{\_}{I}}_{F0} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {\underset{\_}{I}}_{F\_ b} + {\underset{\_}{\; I}}_{F\_ c}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + 0 + 0} \right)} = {\frac{1}{3}{\underset{\_}{I}}_{F\_ a}}}}$

The sequence components are related: I_(F1)=I_(F2)=I_(F0) and finally: I_(F)=I_(F) _(—) _(a)=3I_(F2), thus: a_(F1)=0, a_(F2)=3, a_(F0)=0 (as in the Set I from Table 2)

or I_(F)=I_(F) _(—) _(a)=3I_(F1), thus: a_(F1)=3, a_(F2)=0, a_(F0)=0 (as in the Set II from Table 2)

or I_(F)=I_(F) _(—) _(a)=1,5I_(F1)+1,5I_(F2), thus: a_(F1)=1,5, a_(F2)=1,5, a_(F0)=0 (as in the Set III from Table 2) a-b fault FIGS. 13 a, 13 b: The fault current can be expressed as: I _(F)=I_(F) _(—) _(a) or as:

${\underset{\_}{I}}_{F} = {\frac{1}{2}\left( {{\underset{\_}{I}}_{F\_ a} - {\underset{\_}{I}}_{F\_ b}} \right)}$

Taking into account that in the healthy phase: I_(F) _(—) _(c)=0 and for the faulted phases: I_(F) _(—) _(b)=−I_(F) _(—) _(a) what gives:

${\underset{\_}{I}}_{F1} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {\underset{\_}{a\; I}}_{F\_ b} + {{\underset{\_}{a}}^{2}{\underset{\_}{I}}_{F\_ c}}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {\underset{\_}{a}\left( {- {\underset{\_}{I}}_{F\_ a}} \right)} + {{\underset{\_}{a}}^{2}0}} \right)} = {\frac{1}{3}\left( {1 - \underset{\_}{a}} \right){\underset{\_}{I}}_{F\_ a}}}}$ ${\underset{\_}{I}}_{F2} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {{\underset{\_}{a}}^{2}{\underset{\_}{I}}_{F\_ b}} + {\underset{\_}{a\; I}}_{F\_ c}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {{\underset{\_}{a}}^{2}\left( {- {\underset{\_}{I}}_{F\_ a}} \right)} + {\underset{\_}{a}0}} \right)} = {\frac{1}{3}\left( {1 - {\underset{\_}{a}}^{2}} \right){\underset{\_}{I}}_{F\_ a}}}}$ $\mspace{20mu}{{\underset{\_}{I}}_{F0} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + {\underset{\_}{I}}_{F\_ b} + \; I_{F\_ c}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{F\_ a} + \left( {- {\underset{\_}{I}}_{F\_ a}} \right) + 0} \right)} = 0}}}$

The relation between I_(F1) and I_(F2) is thus:

$\frac{{\underset{\_}{I}}_{F1}}{{\underset{\_}{I}}_{F2}} = {\frac{\frac{1}{3}\left( {1 - \underset{\_}{a}} \right){\underset{\_}{I}}_{F\_ a}}{\frac{1}{3}\left( {1 - {\underset{\_}{a}}^{2}} \right){\underset{\_}{I}}_{F\_ a}} = \frac{\left( {1 - \underset{\_}{a}} \right)}{\left( {1 - {\underset{\_}{a}}^{2}} \right)}}$

Finally:

${\underset{\_}{I}}_{F} = {{\underset{\_}{I}}_{F\_ a} = {{\frac{3}{\left( {1 - {\underset{\_}{a}}^{2}} \right)}{\underset{\_}{I}}_{F2}} = {\left( {1 - \underset{\_}{a}} \right){\underset{\_}{I}}_{F2}}}}$ thus: a _(F1)=0, a _(F2)=1−a, a _(F0)=0 (as in the Set I from Table 2)

or

${\underset{\_}{I}}_{F} = {{\underset{\_}{I}}_{F\_ a} = {{\frac{3}{\left( {1 - \underset{\_}{a}} \right)}{\underset{\_}{I}}_{F1}} = {\left( {1 - {\underset{\_}{a}}^{2}} \right){\underset{\_}{I}}_{F1}}}}$ thus: a _(F1)=1−a , a _(F2)=0, a _(F0)=0 (as in the SET 2 from Table 2)

or

${\underset{\_}{I}}_{F} = {{{0\text{,}5{\underset{\_}{I}}_{F\_ a}} + {0\text{,}5{\underset{\_}{I}}_{F\_ a}}} = {{{\frac{1\text{,}5}{\left( {1 - {\underset{\_}{a}}^{2}} \right)}{\underset{\_}{I}}_{F2}} + {\frac{1\text{,}5}{\left( {1 - \underset{\_}{a}} \right)}{\underset{\_}{I}}_{F1}}} = {{0\text{,}5\left( {1 - \underset{\_}{a}} \right){\underset{\_}{I}}_{F2}} + {0\text{,}5\left( {1 - {\underset{\_}{a}}^{2}} \right){\underset{\_}{I}}_{F1}}}}}$ thus: a _(F1)=0,5(1−a ²), a _(F2)=0,5(1−a), a _(F0)=0 (as in the Set III from Table 2)

(a-b-g) fault, FIG. 14:

${\underset{\_}{I}}_{F} = {{{\underset{\_}{I}}_{F\_ a} - {\underset{\_}{I}}_{F\_ b}} = {{\left( {{\underset{\_}{I}}_{F1} + {\underset{\_}{I}}_{F2} + {\underset{\_}{I}}_{F0}} \right) - \left( {{{\underset{\_}{a}}^{2}{\underset{\_}{I}}_{F1}} + {\underset{\_}{aI}}_{F2} + {\underset{\_}{I}}_{F0}} \right)} = {{\left( {1 - {\underset{\_}{a}}^{2}} \right){\underset{\_}{I}}_{F1}} + {\left( {1 - \underset{\_}{a}} \right){\underset{\_}{I}}_{F2}}}}}$

Thus: a_(F1)=1−a², a_(F2)=1−a, a_(F0)=0 (as in the Sets I, II, III from Table 2)

(a-b-c) or (a-b-c-g) symmetrical faults, FIGS. 15 a, 15 b, 15 c:

Taking the first two phases (a, b) for composing the voltage drop across a fault path one obtains:

${\underset{\_}{I}}_{F} = {{{\underset{\_}{I}}_{F\_ a} - {\underset{\_}{I}}_{F\_ b}} = {{\left( {{\underset{\_}{I}}_{F1} + {\underset{\_}{I}}_{F2} + {\underset{\_}{I}}_{F0}} \right) - \left( {{{\underset{\_}{a}}^{2}{\underset{\_}{I}}_{F1}} + {\underset{\_}{aI}}_{F2} + {\underset{\_}{I}}_{F0}} \right)} = {{\left( {1 - {\underset{\_}{a}}^{2}} \right){\underset{\_}{I}}_{F1}} + {\left( {1 - \underset{\_}{a}} \right){\underset{\_}{I}}_{F2}}}}}$

Thus: a _(F1)=1−a ² , a _(F2)=1−a, a _(F0)=0

Additionally, if a fault is ideally symmetrical the positive sequence is the only component, which is present in the signals. Therefore, we have: a _(F1)=1−a ² , a _(F2)=0, a _(F0)=0 (as in the Sets I, II, III from Table 2). A3. Derivation of the Complex Coefficients in the Fault Current Distribution Factors for the Positive (Negative) Sequence (Table 3) a) Case of the single line with extra link 45 between the substations (FIG. 4)

Let us determine the fault current distribution factor for the positive sequence (the fault current distribution factor for the negative sequence is the same). The equivalent circuit of FIG. 4 with indicated flow of currents for the incremental positive sequence is presented in FIG. 16.

Considering the closed mesh containing: the local segment of the faulted line, the remote segment of the faulted line and the extra link between the substations, one can write down: dZ _(1L) ΔI _(A1)+(1−d)Z _(1L)(ΔI _(A1) −I _(F1))−Z _(1AB) ΔI _(AB1)=0

From the above equation the unknown current from the extra link between the substation can be determined as:

${\Delta\;{\underset{\_}{I}}_{{AB}\; 1}} = {\frac{{\underset{\_}{Z}}_{1L}}{{\underset{\_}{Z}}_{1{AB}}}\left( {{\Delta\;{\underset{\_}{I}}_{A\; 1}} - {\left( {1 - d} \right){\underset{\_}{I}}_{F\; 1}}} \right)}$

Considering the closed mesh containing the source impedances (Z_(1sA), Z_(1sB)) and the extra link (Z_(1AB)) one can write down: Z _(1sA)(ΔI _(A1) +ΔI _(AB1))+Z _(1AB) ΔI _(AB1) +Z _(1sB)(ΔI _(A1) +ΔI _(AB1) −I _(F1))=0

Substituting the previously determined unknown current from the extra link into the above equation one obtains:

${\underset{\_}{k}}_{F\; 1} = {\frac{\Delta\;{\underset{\_}{I}}_{A\; 1}}{{\underset{\_}{I}}_{F\; 1}} = \frac{{{\underset{\_}{K}}_{1}d} + {\underset{\_}{L}}_{1}}{{\underset{\_}{M}}_{1}}}$ where, as in TABLE 3 (Single line, Z_(1AB)≠∝), we have: K ₁ =−Z _(1L) Z _(1AB)−(Z _(1sA) +Z _(1sB))Z _(1L) L ₁ =Z _(1L)(Z _(1sA) +Z _(1sB))+Z _(1AB)(Z _(1L) +Z _(1sB)) M ₁=(Z _(1sA) +Z _(1sB))(Z _(1AB) +Z _(1L))+Z _(1L)Z_(1AB)

If there is no extra link between the substations (Z_(1AB)→∝) one has to consider the closed mesh containing the source impedances (Z_(1sA), Z_(1sB)) and both the segments of the faulted line [dZ_(1L) and (1−d)Z_(1L)]. For this mesh one can write: (Z _(1sA) +dZ _(1L))ΔI _(A1) +[Z _(1sB)+(1−d)Z _(1L)](ΔI _(A1) −I _(F1))=0

After rearrangements one obtains:

${\underset{\_}{k}}_{F\; 1} = {\frac{\Delta\;{\underset{\_}{I}}_{A\; 1}}{{\underset{\_}{I}}_{F\; 1}} = \frac{{{\underset{\_}{K}}_{1}d} + {\underset{\_}{L}}_{1}}{{\underset{\_}{M}}_{1}}}$ where, as in TABLE 3 (Single line, Z_(1AB)→∝), we have: K ₁ =−Z _(1L) L ₁ =Z _(1L) +Z _(1sB) M ₁ =Z _(1sA) +Z _(1sA) +Z _(1L) b) Case of the parallel lines with extra link between the substation (FIG. 5)

Let us determine the fault current distribution factor for the positive sequence (the fault current distribution factor for the negative sequence is the same). The equivalent circuit of parallel lines from FIG. 5 with indicated flow of currents for the incremental positive sequence is presented in FIG. 17.

The healthy parallel line (LB) and the extra link 55 (AB), which are in parallel connection, have been substituted by the equivalent branch with the equivalent impedance:

${\underset{\_}{Z}}_{{{1{LB}}\&}{AB}} = {\frac{{\underset{\_}{Z}}_{1{LB}}{\underset{\_}{Z}}_{1{AB}}}{{\underset{\_}{Z}}_{1{LB}} + {\underset{\_}{Z}}_{1{AB}}}.}$

Considering the closed mesh denoted by (AA, F, BA, BB, AB) one can write: dZ _(1LA) ΔI _(AA1)+(1−d)Z _(1LA)(ΔI _(AA1) −I _(F1))−Z _(1LB&AB) ΔI _(LB&AB1)=0

From the above equation the unknown current from the equivalent branch can be determined as:

${\Delta\;{\underset{\_}{I}}_{{{LB}\;\&}{AB}\; 1}} = {\frac{{\underset{\_}{Z}}_{1{LA}}}{{\underset{\_}{Z}}_{{{1{LB}}\&}{AB}}}\left( {{\Delta\;{\underset{\_}{I}}_{A\;{A1}}} - {\left( {1 - d} \right){\underset{\_}{I}}_{F\; 1}}} \right)}$

Considering the closed mesh containing the source impedances (Z_(1sA), Z_(1sB)) and the equivalent branch (Z_(1LB&AB)) one can write down: Z _(1sA)(ΔI _(AA1) +ΔI _(LB&AB1))+Z _(1LB&AB) ΔI _(LB&AB1) +Z _(1sB)(ΔI _(AA1) +ΔI _(LB&AB1) −I _(F1))=0

Substituting the previously determined unknown current from the healthy line into the above equation one obtains:

${\underset{\_}{k}}_{F\; 1} = {\frac{\Delta\;{\underset{\_}{I}}_{{AA}\; 1}}{{\underset{\_}{I}}_{F\; 1}} = \frac{{{\underset{\_}{K}}_{1}d} + {\underset{\_}{L}}_{1}}{{\underset{\_}{M}}_{1}}}$ where, as in TABLE 3 (Parallel lines), we have: K ₁ =−Z _(1LA)(Z _(1sA) +Z _(1sB) +Z _(1LB&AB)), L ₁ =Z _(1LA)(Z _(1sA) +Z _(1sB) +Z _(1LB&AB))+Z _(1LB&AB) Z _(1sB1) M ₁ =Z _(1LA) Z _(1LB&AB) +Z _(1LA)(Z _(1sA) +Z _(1sB))+Z _(1LB&AB)(Z _(1sA) +Z _(1sB)) where:

${\underset{\_}{Z}}_{{{1{LB}}\&}{AB}} = {\frac{{\underset{\_}{Z}}_{1{LB}}{\underset{\_}{Z}}_{1{AB}}}{{\underset{\_}{Z}}_{1{LB}} + {\underset{\_}{Z}}_{1{AB}}}.}$

In case if the extra link between the substations (Z_(1AB)) is not present one has to substitute: Z_(1LB&AB)=Z_(1LB). 

1. A method to locate a fault from one end of a section of a power line by means of measurements of current, voltage and angles between the phases at a first end of said section, the method comprising: calculating symmetrical components of currents for said current and voltage measurement at said first end, calculating a value of impedance for an extra link between the terminals with the impedance for the positive sequence equal to: $\left( {{\underset{\_}{Z}}_{{{1{LB}}\&}{AB}} = \frac{{\underset{\_}{Z}}_{1{LB}}{\underset{\_}{Z}}_{1{AB}}}{{\underset{\_}{Z}}_{1{LB}} + {\underset{\_}{Z}}_{1{AB}}}} \right)$ where: Z_(1AB)=impedance for the positive sequence of the extra link, Z_(1LA)=positive-sequence impedance of the healthy line, determining a compensation for the shunt capacitance with the aid of an equation of the form: B ₂ ^(comp) ^(—1) (d _(comp) _(—) ₁)² +B ₁ ^(comp) ^(—) ¹ d _(comp) _(—) ₁ +B ₀ ^(comp—) ¹=0 where: B ₂ ^(comp) ^(—) ¹ =A ₂ _(—) _(Re) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ −A ₂ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Re) ^(comp) ^(—) ¹ B ₁ ^(comp) ^(—) ¹ =A ₁ _(—) _(Re) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ −A ₁ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Re) ^(comp) ^(—) ¹ B ₀ ^(comp) ^(—) ¹ =A ₀ _(—) _(Re) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ −A ₀ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Re) ^(comp) ^(—) ¹ determining the zero-sequence current from the healthy line of a section of parallel power lines, calculating a distance to a fault for the parallel line section, calculating the distance to the fault from said first end using a quadratic equation of the form: B ₂ d ² +B ₁ d+B ₀=0 where: B ₂ =A ₂ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₂ _(—) _(Im) A ₀₀ _(—) _(Re) B ₁ =A ₁ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₁ _(—) _(Im) A ₀₀ _(—) _(Re) B ₀ =A ₀ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₀ _(—) _(Im) A ₀₀ _(—) _(Re) and locating the fault utilizing the calculated distances.
 2. The method according to claim 1, wherein the distance to the fault is calculated using an equation of the form: $\begin{matrix} {{{{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{1L}d^{2}} + {\left( {{{\underset{\_}{L}}_{1}{\underset{\_}{Z}}_{1L}} - {{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{AA\_ p}}} \right)d} - {{\underset{\_}{L}}_{1}{\underset{\_}{Z}}_{AA\_ p}} + {R_{F}{\underset{\_}{M}}_{1}\frac{\left( {{{\underset{\_}{a}}_{F\; 1}\Delta\;{\underset{\_}{I}}_{{AA}\; 1}} + {{\underset{\_}{a}}_{F\; 2}{\underset{\_}{I}}_{{AA}\; 2}}} \right)}{{\underset{\_}{I}}_{AA\_ p}}}} = 0} & (8) \end{matrix}$ where: ${\underset{\_}{Z}}_{AA\_ p} = \frac{{\underset{\_}{V}}_{AA\_ p}}{{\underset{\_}{I}}_{AA\_ p}}$ calculated fault loop impedance.
 3. The method according to claim 1, wherein the distance to the fault is calculated using an equation of the form: A ₂ d ² +A ₁ d+A ₀ +A ₀₀ R _(F)=0 where: A ₂ =A ₂ _(—) _(Re) +jA ₂ _(—) _(Im) =K ₁ Z _(1LA) A ₁ =A ₁ _(—) _(Re) +jA ₁ _(—) _(Im) =L ₁ Z _(1LA) −K ₁ Z _(AA) _(—) _(p) A ₀ =A ₀ _(—) _(Re) +jA ₀ _(—) _(Im) =−L ₁ Z _(AA) _(—) _(p) ${A_{00{\_ Re}} + {j\; A_{00{\_ Im}}}} = \frac{{\underset{\_}{M}}_{1}\left( {{{\underset{\_}{a}}_{F1}\Delta\;{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{AA2}}} \right)}{{\underset{\_}{I}}_{AA\_ p}}$ ${\underset{\_}{Z}}_{AA\_ p} = {\frac{{\underset{\_}{V}}_{AA\_ p}}{{\underset{\_}{I}}_{AA\_ p}} = {{calculated}\mspace{14mu}{fault}\mspace{14mu}{loop}\mspace{14mu}{impedance}}}$ K₁, L₁, M₁=coefficients gathered in TABLE
 3. 4. The method according to claim 1, further comprising: determining source impedance at said first end as a representative value, and determining a value for source impedance at said second end as a representative value.
 5. The method according claim 1, further comprising calculating symmetrical components of currents for said current and voltage measured at said first end by: inputting instantaneous phase voltages, filtering the values to determine the phasors, and calculating phasors of symmetrical components of voltages.
 6. The method according to claim 1, further comprising calculating symmetrical components of currents for said current and voltage measured at said first end by: inputting instantaneous phase currents and instantaneous zero-sequence current from a healthy line, filtering the values to determine the phasors, and calculating phasors of symmetrical components of currents.
 7. The method according to claim 1, further comprising determining a compensation for shunt capacitance by means of an equation of the form: ${{{\underset{\_}{A}}_{2}^{{comp\_}1}\left( d_{{comp\_}1} \right)}^{2} + {{\underset{\_}{A}}_{1}^{{comp\_}1}d_{{comp\_}1}} + {\underset{\_}{A}}_{0}^{{comp\_}1} + {{\underset{\_}{A}}_{00}^{{comp\_}1}R_{F}}} = 0$ ${\underset{\_}{A}}_{2}^{{comp\_}1} = {{A_{2{\_ Re}}^{{comp\_}1} + {j\; A_{2{\_ Im}}^{{comp\_}1}}} = {{{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{1L}^{long}{\underset{\_}{A}}_{1}^{{comp\_}1}} = {{A_{1{\_ Re}}^{{comp\_}1} + {j\; A_{1{\_ Im}}^{{comp\_}1}}} = {{{{\underset{\_}{L}}_{1}{\underset{\_}{Z}}_{1L}^{long}} - {{\underset{\_}{K}}_{1}{\underset{\_}{Z}}_{A\_ p}^{{comp\_}1}{\underset{\_}{A}}_{0}^{{comp\_}1}}} = {{A_{0{\_ Re}}^{{comp\_}1} + {j\; A_{0{\_ Im}}^{{comp\_}1}}} = {{{- {\underset{\_}{L}}_{1}}{\underset{\_}{Z}}_{A\_ p}^{{comp\_}1}{\underset{\_}{A}}_{00}^{{comp\_}1}} = {{A_{00{\_ Re}}^{{comp\_}1} + {j\; A_{00{\_ Im}}^{{comp\_}1}}} = \frac{{\underset{\_}{M}}_{1}\left( {{{\underset{\_}{a}}_{F1}\Delta\;{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{AA2}}} \right)}{{\underset{\_}{I}}_{A\_ p}^{{comp\_}1}}}}}}}}}$ (21a) where: ${\underset{\_}{Z}}_{A\_ p}^{{com\_}1} = \frac{{\underset{\_}{V}}_{A\_ p}}{{\underset{\_}{I}}_{A\_ p}^{{comp\_}1}}$ fault loop impedance calculated from: V _(A) _(—) _(p)—original (uncompensated) fault loop voltage, I_(A) _(—) _(p) ^(comp) ^(—) ¹=a₁I_(A1) _(—) _(comp) _(—) ₁+a₂I_(A2) _(—) _(comp) _(—) ₁+a₀I_(A0) _(—) _(comp) _(—) ₁—fault loop current composed of the positive, negative and zero sequence currents obtained after deducing the respective capacitive currents from the original currents, and K₁, L₁, M₁=coefficients gathered in TABLE
 3. 8. The method according to claim 1, further comprising measuring the source impedance Z_(1sA) at said first end A.
 9. The method according to claim 1, further comprising: measuring the source impedance Z_(1sB) at said second end, sending a communication of the measured value of source impedance Z_(1sB) at said second end to a fault locator at said first end.
 10. The method according to claim 1, further comprising determining the distance to a single phase to ground fault without measurements from an operating healthy parallel line by means of complex coefficients P₀ according to a formula of the form: ${\underset{\_}{P}}_{0} = \frac{{\underset{\_}{Z}}_{0{LB}} - {\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{0{LA}} - {\underset{\_}{Z}}_{0m}}$ and K₁, L₁, M₁ according to K ₁ =−Z _(1LA)(Z _(1sA) +Z _(1sB) +Z _(1LB)) L ₁ =−K ₁ +Z _(1LB) Z _(1sB) M ₁ =Z _(1LA) Z _(1LB) +Z _(1LA)(Z_(1sA) +Z _(1sB))+Z _(1LB)(Z _(1sA) +Z _(1sB)).
 11. The method according to claim 1, further comprising determining the distance to a single phase to ground fault without measurements from switched off and grounded parallel line by means of complex coefficients P₀ according to ${\underset{\_}{P}}_{0} = {- \frac{{\underset{\_}{Z}}_{0{LB}}}{{\underset{\_}{Z}}_{0m}}}$ and K₁, L₁, M₁ according to K ₁ =−Z _(1LA) L ₁ =Z _(1LA) +Z _(1sB) M ₁ =Z _(1sA) +Z _(1sA) +Z _(1LA).
 12. The method according to claim 1, further comprising determining the distance to a single ground fault using a first order formula of the form: $d = {\frac{{imag}\left\{ {{\underset{\_}{V}}_{AA\_ p}\left\lbrack {3\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}{{imag}\left\{ {\left( {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ p}} \right)\left\lbrack {3\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}.}$
 13. The method according to claim 1, further comprising determining the distance to a phase-to-phase ground fault using pre-fault measurements and a first order formula of the form: $d = {\frac{{imag}\left\{ {{\underset{\_}{V}}_{AA\_ p}\left\lbrack {\underset{\_}{W}\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}{{imag}\left\{ {\left( {{\underset{\_}{Z}}_{1{LA}}{\underset{\_}{I}}_{AA\_ p}} \right)\left\lbrack {\underset{\_}{W}\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)} \right\rbrack}^{*} \right\}}.}$
 14. The method according to claim 1, further comprising determining the distance to a phase-to-phase ground fault avoiding pre-fault measurements and using a first order formula of the form: $d = {\frac{{imag}\left\lbrack {\left( {{\underset{\_}{V}}_{a} + {\underset{\_}{V}}_{b}} \right)\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)^{*}} \right\rbrack}{{imag}\left\lbrack {{{\underset{\_}{Z}}_{1{LA}}\left( {{\underset{\_}{I}}_{a} + {\underset{\_}{I}}_{b} + {2{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{{AA}\; 0}2{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{{AB}\; 0}}} \right)}\left( {{\underset{\_}{I}}_{{AA}\; 0} - {{\underset{\_}{P}}_{0}{\underset{\_}{I}}_{{AB}\; 0}}} \right)^{*}} \right\rbrack}.}$
 15. A device for locating a fault from one end of a section of a power line having means for receiving and storing measurements of current, voltage and angles between the phases at one first end, means for receiving and storing a detection of a fault condition between said first and second ends, the device comprising: means for calculating symmetrical components of currents for said current and voltage measured at said first end, means for calculating a value of impedance for an extra link between the terminals, means for determining a compensation for shunt capacitance, means for determining the zero-sequence current from the healthy line of a section of parallel power lines, means for calculating a distance to a fault for the parallel line section, means for calculating a distance from said first end to the fault.
 16. The device according to claim 15, further comprising: means for determining a value for source impedance at said first end, means for determining a value for source impedance at said second end.
 17. The device according to claim 15, further comprising: means for receiving a measurement of source impedance at said first end.
 18. The device according to claim 15, further comprising: means for receiving a measurement of source impedance made at said second end.
 19. The device according to claim 15, further comprising means to receive a measured value for remote source impedance at said second end communicated by means of a communication channel.
 20. Use of a fault locator device according to claim 15, by a human operator to supervise a function in an electrical power system.
 21. Use of a fault locator device according to claim 15, by means of a process running on one or more computers to supervise and/or control a function in an electrical power system.
 22. Use of a fault locator device according to claim 15, to locate a distance to a fault in a power transmission or distribution system.
 23. Use of a device according to claim 15, for locating a fault on parallel power lines.
 24. A computer program product, comprising a computer readable medium; and computer code means and/or software code portions recorded on the computer readable medium for making a computer or processor perform a method for locating a fault from one end of a section of a power line by means of measurements of current, voltage and angles between the phases at a first end of said section, the method comprising calculating symmetrical components of currents for said current and voltage measurement at said first end, calculating a value of impedance for an extra link between the terminals with the impedance for the positive sequence equal to: $\left( {{\underset{\_}{Z}}_{{{1\;{LB}}\&}{AB}} = \frac{{\underset{\_}{Z}}_{1{LB}}{\underset{\_}{Z}}_{1{AB}}}{{\underset{\_}{Z}}_{1{LB}} + {\underset{\_}{Z}}_{1{AB}}}} \right)\mspace{14mu}{where}\text{:}$ Z_(1AB)=impedance for the positive sequence of the extra link, Z_(1LA)=positive-sequence impedance of the healthy line, determining a compensation for the shunt capacitance with the aid of an equation of the form: B ₂ ^(comp) ^(—) ¹(d _(comp) _(—) ₁)² +B ₁ ^(comp) ^(—) ¹ d _(comp) _(—) ₁ +B ₀ ^(comp) ^(—) ¹=0 where: B ₂ ^(comp) ^(—) ¹ =A ₂ _(—) _(Re) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ −A ₂ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(Re) ^(comp) ^(—) ¹ B ₁ ^(comp) ^(—) ¹ =A ₁ _(—) _(Re) ^(comp) ^(—) ¹ A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ −A ₁ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(Re) ^(comp) ^(—) ¹ B ₀ ^(comp) ^(—) ¹ =A ₀ _(—) _(Re) ^(comp) ^(—) ¹ −A ₀₀ _(—) _(Im) ^(comp) ^(—) ¹ A ₀ _(—) _(Im) ^(comp) ^(—) ¹ A ₀₀ _(Re) ^(comp) ^(—) ¹ determining the zero-sequence current from the healthy line of a section of parallel power lines, calculating a distance to a fault for the parallel line section, calculating the distance to the fault from said first end using a quadratic equation of the form: B ₂ d ² +B ₁ d+B ₀=0 where: B ₂ =A ₂ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₂ _(Im) A ₀₀ _(—) _(Re) B ₁ =A ₁ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₁ _(Im) A ₀₀ _(—) _(Re) B ₀ =A ₀ _(—) _(Re) A ₀₀ _(—) _(Im) −A ₀ _(Im) A ₀₀ _(—) _(Re) and locating the fault utilizing the calculated distances. 